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The Sliding Surface* 
ByreG.. la bINCH 


34th Guthrie Lecture, delivered 3rd March 1950 


$1. INTRODUCTION 


IFE 1s possible because some assemblages of matter can store sunlight as 
is chemical energy. Other organisms, unable to do this directly, absorb these 
stores to supply their own energy needs; and some, like man, prey on both 
to live. Man alone supplements the body’s production of physical energy by 
using inanimate devices to unlock the stores of energy latent in the remnants of 
past life. The key is the machine, which works because smooth surfaces can 
slide over each other. The sliding surface is, therefore, fundamental to our 
modern way of life. 

The sliding surface was known long before the machine; away back in 
antiquity man had used it to ease the burden on his back. He had also learnt 
the value of a generous supply of some slippery substance. So when the machine 
did arrive the engineer saw to it that its sliding surfaces were lubricated; so well, 
indeed, that they were kept apart by an oil layer many hundreds or even thousands 
of molecules thick. Nowadays, however, sliding surfaces are usually so heavily 
pressed together that the lubricant film cannot be more than a few molecules 
thick, and is often penetrated or disrupted so that the sliding surfaces touch in 
places. 

Thick-layer, or, as it is generally known, hydrodynamic lubrication (Reynolds 
1886) favours easy sliding without wear other than that due to corrosion and erosion, 
but with boundary (or molecular film) conditions + of lubrication, where the sliding 
surfaces often come into contact, suppression of wear becomes the main problem. 
Thus, in hydrodynamic lubrication chief interest is directed to the lubricant, but 
in boundary lubrication it is the sliding surface which claims first attention. 
What causes rubbing surfaces to wear and how do they wear? What peculiarity 
in the surface helps to reduce wear and how does it work? How do the surfaces 
slide when they touch each other? What is it about the surfaces that helps them 
to slide and how does it help? We must first seek the answers to questions like 
these before turning to the study of the lubricant boundary film. 

As far back as the eighteenth century Coulomb (1781) thought that resistance 
to sliding was due to the interlocking of surface asperities, and that the fracture of 
these during sliding was responsible for wear. Ewing (1912) and Hardy and 

* This Lecture will also be published in Section A in August. 
+ The term ‘ boundary conditions’ was first used by Reynolds. 
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Hardy (1919), on the other hand, believed that the interplay of surface cohesive 
forces was the dominant factor. While Coulomb’s view may once have served 
a useful purpose, it is completely at variance with modern engineering experience. 
Furthermore, the results of critical surface studies by Beilby (1921), Hardy (1936), 
Finch (1937), Bowden and Tabor (1939) and Parker and Hatch (1950) do not 
support Coulomb’s hypothesis but are readily interpreted in terms of surface 
cohesive forces. 

The forces which produce cohesion within a block of metal must also act at the 
surface of contact between two blocks; they should, therefore, resist stress as if 
the blocks were one. There are, indeed, facts which support this conclusion. Slip 
gauges left ‘wrung’ together eventually seize and are damaged when separated. 
Hardy and Hardy (1919) found that some solids seized together on contact. 
Fifteen years ago I showed—and it must have been known long before then— 
that mica cleavage surfaces could adhere firmly when brought together (Demonstra- 
tion). But experience also shows that surfaces can, and often do, slide easily 
over each other, even in the absence of lubricating oil. 

Although a solid exerts a cohesive force on an atom beyond its surface, the 
intensity of the force decreases rapidly with distance, to become negligible a few 
atomic radii away. Contact between two surfaces, therefore, implies approach 
to within this distance. It is possible, by cleaving certain mineral crystals, to 
produce atomically smooth plane surfaces of considerable area, but otherwise 
such a degree of smoothness cannot in general be achieved or even approached. 
Hence, when two nominally plane surfaces formed, let us suppose, on two blocks 
of the same metal are placed together, contact can occur only at three points. 
If the surfaces are now loaded, stress at the points of contact must exceed the 
elastic limit. Plastic flow will then enlarge the points of contact to areas, and fresh 
points of contact expanding to still more areas will be established, until the load 
is borne without further deformation of the material.* A considerable effort 
should now be required to pull or shear the blocks apart because, at the areas of 
contact, they should have become united. Experiment, however, gives a different 
result. Thus, mica cleavage surfaces which have been exposed to the atmosphere 
for some time, or momentarily immersed in water, will no longer adhere, no 
matter how hard they are forced together (Demonstration). If, instead of being 
‘wrung’, two slip gauges are merely pressed together, they do not stick but fall 
apart when the pressure is released (Demonstration). Mechanically or anodically 
polished surfaces of copper, iron, silver, tin and other metals behave similarly, 
provided the pressure does not cause general plastic flow at the surfaces 
(Demonstration). An atomic degree of smoothness, superimposed, it is true, 
on a large-scale waviness of small amplitude, may be imparted to a face of a 
copper single-crystal by anodic polishing; nevertheless, a degree of pressure 
short of that causing permanent distortion will not even partially weld two such 
surfaces together (Demonstration). In this case the undimmed clarity of the 
electron diffraction pattern shows that any adsorbed layer cannot be more than two 
or three atoms thick (Figure 1).+ 

* I frame these conclusions on the basis of Amontons’ (1699) experimental law according to 
which the coefficient of friction is independent of the area of the sliding surfaces. The validity of 
this law has often been confirmed. ‘The reasons for apparent exceptions are well understood ; they 


do not affect the issue. 
} Figures 1 and 3-20 are photographs and are printed together opposite p. 472. 


34th Guthrie Lecture 467 


It may therefore be inferred that true contact is not as a rule established when 
even atomically smooth solid surfaces are pressed together, but that they are 
kept apart by adsorbed surface layers which, while satisfying the cohesive forces. 
at the surfaces, exert little or no attraction towards each other. The remarkable 
resistance of these layers to penetration appears to be independent of their nature 
(if it is assumed that this varies with the metal) and determined mainly by the 
physical properties of the metal, general plastic flow of which might be expected 
to rupture such a layer by extending its surface, and thus allowing welding to. 
occur through direct metal-to-metal contact. It will be evident that this adsorbed 
layer is of vital importance. Without it, a pair of surfaces at rest under load 
would seize together over an area proportional to the imposed load. It cannot be 
supposed that welds are formed only to be broken when, on release of pressure,. 
the surfaces are forced apart by a recovery of the metal about those areas where 
the elastic limit has not been exceeded; for the metal in such welds, being itself 
compressed to the elastic limit, would have the greater range of recovery on release 
of the pressure. 

When the surfaces are not only pressed together but forced to slide past each 
other, the effectiveness of the adsorbed layers in preventing metal-to-metal 
contact varies greatly according to the nature of the metal. Anodically polished 
single-crystal copper or aluminium surfaces at first slide then suddenly lock 
firmly together, even with light loading (Demonstration). Similarly prepared 
polycrystalline surfaces of these metals also seize quickly, though not quite so 
readily as single crystals (Demonstration); the corresponding mechanically 
polished surfaces slide more easily and are much more resistant to seizure 
(Demonstration). Anodically polished polycrystalline silver or iron single- 
crystal pairs slide freely and must be heavily loaded before signs of seizure appear 
(Demonstration), and similar experiments with combinations of anodically 
polished single-crystal and polycrystalline surfaces of various metals reveal 
further wide differences in the ease or otherwise with which sliding and particularly 
seizure take place. But whether this is to be attributed to the physical properties 
of the metal, as in the case of the loaded surface pair at rest, or whether the nature 
of the adsorbed layer is also significant, remains to be seen. 

The existence of these layers has long been known. Faraday (1838) and 
De la Rive (1839) recognized their importance in heterogeneous catalysis, and they 
have since been much studied by the chemist. ‘The physicist’s interest seems to: 
have been first fully aroused by their effect on thermionic and photoelectric 
emission (Langmuir 1916). The engineer, so long accustomed to hydrodynamic 
lubrication, has been slower to appreciate their value. 

The layer is difficult to remove from a metal surface. It can be done by 
scraping or evaporating im vacuo and so exposing a fresh surface, or by heating 
to a high temperature 7m vacuo or in the presence of atomic hydrogen ; but the layer 
re-forms spontaneously and rapidly, even at room temperature, on exposure to air. 
These facts suggest that the adsorption is chemical rather than physical, and that 
the layer is an oxide of the metal. This has been amply confirmed by electron 
diffraction. 

Although oxidation is the immediate result of the exposure of a fresh metal 
surface to the atmosphere, the microbalance reveals that a moisture film of consider- 
able thickness (as much as 80 molecular layers and possibly more) may be adsorbed 
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on the oxide surface. When fresh, this layer is as a rule not strongly held and is 
removed by the electron beam in the evacuated diffraction camera. 


§2. THE SLIDING ELEMENT 

It is clear that the sliding element of the modern machine must be a highly 
complex system. Actually its simplest form comprises at least nine distinct 
components (Figure 2). These are (i) and (ii) the metallic backings supporting 
(iii) and (iv) the pair of metal surface layers or so-called bearing surfaces which, in 
turn, carry (v) and (vi) their appropriate oxide layers each with its (vil) and (vii1) 
adsorbed boundary layer of lubricant, and (ix) an oil film of varying thickness 
held between the boundary layers and more or less charged with debris, moisture 
and gas bubbles. The components in each pair may differ widely in nature and 
properties, even in the case of the two lubricant boundary layers. Furthermore, 
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Figure 2. Diagram of the sliding element. 
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a change in any one component may also affect the behaviour of the others. It is 
this astonishing complexity in what at first sight appears to be a simple system that 
explains why, in the past, research has been for the most part directed towards the 
special circumstance of practice, rather than to a fundamental study from which 
broad generalizations might be drawn. Clearly, the physicist’s approach to the 
problem must be to break down the sliding element into its simplest components 
and to examine each in turn, first by itself and then in relation to the others. 


§3. THE BACKING 

The backing is rigidly secured, directly or indirectly, to an appropriate member 
of the machine. Its function is to support the bearing surface so that its motion 
is constrained to the desired degree or degrees of freedom with respect to the other 
sliding member. In the case of a journal bearing, the force tending to rotate the 
backing in its housing may be considerable and is given by the product of the 
coefficient of friction between the sliding members and the load. Since the polar 
distribution of the load may vary periodically between wide extremes, particularly 
in the case of the internal combustion engine, the backing tends to oscillate in the 
housing. ‘This may lead to serious frettage corrosion and loss of fit. 

Backing and bearing surface may be composed of the same metal, as in the 
shaft of a journal bearing; or the metals may differ, as in the liner of the housing 
within which the shaft rotates. As the properties of the backing may affect the 
behaviour of the bearing surface, it will be convenient to defer its further 
discussion. 
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§4. THE BEARING SURFACE 


‘The bearing surface is developed in practice as a plane or surface of revolution 
by abrasion processes such as cutting, grinding, honing and lapping. These 
operations reduce the crystal size of the metal at the surface, and often toa consider- 
able depth below it (Figure 3). This ‘work-hardening’ can seriously impair 
the efficiency of the bearing surface, and for this reason the more severe operations 
of cutting and grinding should be followed by annealing to restore the uniformity 
of crystal size and relieve stress (Finch 1945). The texture of the finished surface 
is, in terms of atomic dimensions, exceedingly rough and made up of rows of steep 
hills and valleys, each with smaller irregularities and lying roughly parallel to 
the direction of cutting. The ridge crests are seen by electron diffraction to be 
crystalline (Figure 4). 

General experience has shown that the fresh bearing surface must not be 
subjected to full loads and sliding speeds until it has been ‘run-in’ under light 
loads with copious supplies of lubricant. If the protuberances are not gradually 
reduced, as they are by ‘running-in’, the crystal structure will be broken up to such 
a depth below the surface as permanently to impair the sliding and resistance to 
wear of the surface. The effect of ‘running-in’ may be explained briefly as 
follows. Consider, first, bearing surfaces of similar metal. When, through 
breakdown of the oil and oxide layers, metal-to-metal contact is established between 
opposing high spots of the bearing surfaces, a weld is formed. As sliding motion 
proceeds the weld is stressed and eventually ruptured. If the crystal structure in 
the surfaces has not previously been broken up and therefore hardened, the weld 
will shear at or near one of the surfaces. On the other hand, if, owing to previous 
unduly severe conditions of machining or sliding, work-hardening has occurred to 
some depth below the surface, the weld, instead of rupturing virtually in the plane 
of sliding, will continue to adhere to one surface and pull out work-hardened metal 
to a considerable depth from the other (Figure 5). The protruding metal will 
then score and break up the crystal structure in this surface until it has been worn 
away; or it may break off and roll between the surfaces and damage them by 
crushing (Figure 6). 

It is now known that sliding causes a profound and important change, not only 
in the texture but also in the structure of the metal surface. At first, rubbing of 
the sharper asperities on the opposing fresh bearing surfaces results largely in 
wear and the production of debris. As the surface becomes smoother, however, 
wear decreases and much of the metal removed from the hills ‘flows’, as Beilby 
(1921) puts it, into the valleys. ‘Thus ‘running-in’ levels off the mountain ridges 
and fills the valleys to leave a texture like undulating hills. ‘This smoothing is 
accompanied by a remarkable change in structure (Finch, Quarrell and Wilman 
1935). Not only are the crystals broken up to a depth below the surface depending 
on the severity of ‘running-in’, but the surface layer structure becomes 
amorphous (Figure 7), as if ‘flowing’ consisted of the transfer of molten metal 
from one part of the surface to another where, having escaped the polishing 
action, it freezes too rapidly to crystallize. ‘The thickness of this so-called Beilby 
layer depends, like that of the work-hardened crystalline zone beneath, on the 
severity of the ‘running-in’ conditions. Layers of only about fifteen atoms 
thick have been produced (Hopkins 1935) and others about 100 atoms thick 
have been successfully detached from their substrates (Cochrane 1938); but a 
vigorous polishing results in a much thicker layer (Finch, Quarrell and Roebuck 
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The following experiments will show how great is the protection afforded to 
the oxide layer by the support of even a very thin Beilby layer. Severe welding 
and extensive damage occurs when two anodically polished copper single-crystal 
surfaces are slid past each other (Figure 8) (Demonstration). The behaviour of a 
similar pair of surfaces, on one of which, however, a very thin Beilby layer has 
been formed by light polishing with rouge, is very different. Much heavier 
pressure is needed to cause welding, and on forcible separation of the surfaces 
metal is pulled out from the polished surface only to a shallow depth (Figure 9) 
(Demonstration). Thus the Beilby layer performs important functions in 
sliding. Its comparative smoothness favours a wide distribution of load where 
‘high spots’ meet as the bearing surfaces slide past each other; it is more rigid than 
the underlying crystalline metal and thus better able to support the oxide layer and 
preserve it from injury; and a cold weld, with its tendency to pull out near the 
zone where the amorphous merges into the undisturbed crystalline structure, will 
do least damage if the Beilby layer is thin. 

The relative amounts of metal worn away (and thus lost to the surface) and 
metal which is merely displaced to form the Beilby layer, are affected by the 
direction of the ridge and valley lines produced in machining in relation to that of 
sliding (Finch and Zahoorbux 1937). It has been found that the heaviest wear 
occurs when these lines in both bearing surfaces lie in the direction of sliding, and 
that wear is least when they are at right angles to the sliding direction. Ifthe lines 
on one surface are parallel to the direction of sliding, while those on the other are 
at right angles to it, then wear of the former surface is high. The explanation is 
simple. When the lines in both surfaces are parallel to the sliding direction, the 
metal, severely stressed in localized contacts has farther to flow before it can escape 
and congeal beyond the pressure zone, and is thus more likely to form unattached 
detritus than when the ridges and valleys lie at right angles to the sliding. The 
importance of the ‘run’ of tooling tracks in relation to the sliding direction is now 
appreciated in engineering practice. 

We have seen that when two nominally plane surfaces are pressed together 
the total cross-sectional area of the resulting intermetallic contact is proportional 
to the load. The force required to start sliding is therefore also proportional 
to the load; its value per unit load is known as the coefficient of static friction. 
To maintain sliding, work must be done. ‘The necessary force is again propor- 
tional to the load, provided shear occurs at the interface and the strength of the 
metal in the contacts is unchanged. Under such conditions the dynamic and 
static coefficients of friction should be equal. Micrographic examination shows, 
however, that the cold welds formed on contact between, for example, single- 
crystal copper surfaces are strain-hardened as sliding proceeds, and also that the 
surrounding metal is similarly affected to a considerable depth (Figures 5 and 8). 
Hence, instead of the welds shearing at the interface, metal is pulled out bodily 
from below one surface to be dragged along and ruck and score the other (Figure 5). 
Under such conditions the dynamic coefficient exceeds the static. In normal 
practical conditions of sliding, of course, the reverse is the case, as the oxide and 
lubricant films together greatly limit intermetallic contact during sliding and may 
even in some cases inhibit it. Nevertheless, as a rule, spasmodic intermetallic 
contact does occur, and its contribution to the total resistance to sliding may be 
called intermetallic friction. Its most undesirable feature is that it causes wear. 
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A sliding element in which the intermetallic component of the total friction is 
considerable is therefore unsatisfactory. 

So far we have considered bearing surface pairs of the same metal. In practice, 
however, where sliding speeds are high, this combination is avoided except in 
the case of cast iron. Otherwise the engineer’s empirical rule is that the mating 
bearing surfaces should consist of metals of widely differing melting points. 
It can, indeed, be inferred from Bowden and Ridler’s (1936) measurements of the 
temperatures attained at metal contacts in sliding that localized fusion followed by 
hot welding should occur more readily with like than with unlike metals. The 
matter, however, goes deeper than this. Solid mercury and gold would make an 
unsatisfactory pair of bearing surfaces. 

Now, we have seen that in addition to hot welding due to melting of rubbing 
high spots on both bearing surfaces, cold welding also occurs on contact of the solid 
metals. Indeed, under the conditions of the experiments I showed you, the 
friction and surface damage was entirely due to cold welding and its after effects. 
Even with like metals, it is only under severe sliding conditions and more particu- 
larly at high speeds that hot welding becomes important. It will be evident, 
therefore, that in selecting metals for a bearing surface pair, one criterion is the 
intensity of the cohesive force across the interface of intermetallic contact, rather 
than disparity in their melting points. If this interfacial cohesion is less than that 
within the weaker of the two metals, intermetallic contacts will shear in the inter- 
face. ‘The other criterion, which concerns the nature of the oxide layer, will be 
discussed later. 

The problem of estimating this cohesive force from a knowledge of the 
structures of the outer electron shells of the atoms of which the contacting surfaces 
are composed has not been solved. However, a study of the coordination relation- 
ships between epitaxial metallic deposits and metallic single-crystal substrates is 
now being carried out in my laboratory (Finch and Sun 1936, Finch and Williams 
1937, Finch, Wilman and Yang 1947); this seems to offer a promising line of 
experimental approach to the problem. In the meantime, as a rough and ready 
rule, we may take it that the more rigid the metals of a pair and the higher their 
melting points, the less will be the tendency to form cold welds.* 

It would, I think, be fair to say that, as in the development of most aspects 
of the sliding surface, the choice of metals for bearing surface pairs has been 
empirical. The modern journal bearing practice is to use pairs in which one metal 
is hard and has a high melting point and the other is relatively soft and fusible. 
In the sliding of such combinations new phenomena appear. 

Bowden and Hughes (1937) showed that when certain solids were rubbed 
together it was not the harder surface which caused the softer to flow and thus be 
polished, but that the more fusible surface was polished by the more refractory 
one, irrespective of their relative hardnesses at room temperature. In practice, 
however, sliding surfaces may well behave differently. For example, a type of 
journal bearing in common use consists of a superficially hardened alloy steel 


* It would be virtually impossible to put together even two perfectly plane complementary single- 
erystal surfaces on a pair of blocks of the same metal in such a way as to restore the original atom 
coordination across the interface. The higher the mobility of the atoms, however, the greater is the 
chance of the atoms in the opposing faces moving into and being retained in sites of lowest potential 
energy. This is akin to a process of recrystallization whereby the two blocks eventually become 
welded together into a single crystal. 
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shaft rotating against a bearing surface of a rather soft tin-base alloy melting at 
about 250°c. Since the steel melts at roughly 1,400° c. the alloy surface becomes — 
polished as one expects; what is surprising, however, is that the hardened surface 

of the shaft also flows and is polished as running-in of the bearing proceeds. 
Similar observations have been made in reciprocating sliding. Indeed, in all 
practical cases so far examined this mutual polishing of sliding metal surfaces has _ 
proved to be the rule, irrespective of the difference in hardness or fusibility of the 
bearing surface metals. The explanation is to be sought in the nature and 
properties, including the geometry, of the oxide layers and metals of the bearing 
surfaces, and will be discussed later. 

Quite apart from the reduction in welding, and hence in friction and wear, 
further advantages can accrue from the use of a relatively fusible soft alloy sliding 
on a hard refractory bearing surface. Owing to the comparative ease with which — 
the softer metal is plastically deformed, it was long thought necessary to restrict 
the load intensity. The further disadvantage of the low ductility and proneness 
to fracture of the customary rather thick strip of soft metal was overcome, at least 
in part, by welding to a massive backing of a metal of more suitable mechanical 
properties, which was then secured in an appropriate housing. ‘The resulting 
bearing, however, was large and massive in relation to the applied load. Hopkins 
and Palm (1932) realized that, by taking advantage of the cohesion between the 
bearing surface and the backing metals, the softer metal could be made to with- 
stand a load intensity far greater than hitherto thought practicable. ‘The under- 
lying principle is simple. The ability of a soft body to resist deformation when 
compressed between two hard plane surfaces increases as its thickness is reduced 
by plastic flow until, long before the soft material can be totally extruded, the 
hard metals themselves begin to flow and deform. For example, it is impos- 
sible to force out completely even an oil film from between two plane surfaces, 
no matter how heavy the load or how rigid the metal may be. ‘This is due to 
the reinforcing of the cohesion between the oil molecules by their adhesion to 
the surfaces. In fact, the thinner the film, the more do its mechanical properties 
approximate to those of the solid substrate. Similarly, whilst a bearing metal 
layer of the once customary thickness can support only a relatively low load 
intensity, a much thinner layer can carry a far heavier load without being 
permanently deformed. ‘The thickness of this layer need be no greater than that 
corresponding to the maximum wear loss which can be tolerated before the fit 
becomes so slack as to render the bearing useless. For example, the ‘ big-end’ 
crankshaft bearing of an internal combustion engine, which previously consisted 
of a 5/100in. thick tin-base alloy fused on to a massive bronze backing, has now © 
been replaced by a 5/1,000 in. thick layer on a thin steel strip, the weights of the 
liners being 218 and 60 gm. respectively (Exhibits). 

The tin-base type of bearing metal alloy, though relatively soft, is too hard to 
be scratched with the finger nail. ‘This form of alloy, though extremely wear- 
resisting, causes considerable wear of the hardened steel surface against which it 
slides. As we shall see later, this is due partly to the softness of the alloy and partly 
to the nature of its oxide layer. ‘There is little evidence to support the view that 
embedded detritus particles are the cause. A further notable advance therefore 
was made by the recent introduction of bearing metals whose surfaces cause very 
little wear of the steel bearing surface. Lead and indium, the chief metals of this 
kind, are so soft that they and their alloys can be scratched with the finger nail. 
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‘igure 1. Electron-diffraction pattern from a copper 
single-crystal electropolished cube face; beam along 
cube edge. The sharp Kikuchi lines and the narrow 
circular Laue zone show that the surface is atomi- 
cally smooth and of quite remarkable lattice per- 
fection. The clarity in detail near the shadow edge 
is such as to show that the oxide layer cannot be 
more than two or three molecules thick. 


Figure 3. An iron single-crystal (110) face was 
lightly honed along a [110] direction and then 
etched to expose the sub-surface damage. 
The pattern (beam along [001]) reveals that 
the original single-crystal structure has been 
distorted by a process of rotational slip 
(Wilman 1950) on (001) about a [001] axis. 
This slipping represents the initial stage in 
crystal break-up. Note that the fragments 
are still coherent across the planes of rota- 
tional slip. 
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Figure 4. -«-iron pattern from the 
freshly honed surface of an austen- 
itic steel cylinder liner. 


Figure 5. Scoring damage due to welding between a 
pair of copper single-crystal sliding surfaces. 
450. 
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Figure 17. The same surface after 30 hours’; 
exposure to air at room _ temperature| 
becomes covered with an amorphous; 
oxide layer which is thick enough (>40a.)| 
to obscure the single-crystal pattern. 


Figure 16. A copper single-crystal surface immediately 
after anodic polishing. 


Figure 18. Cuprous oxide growing Figure 19. y—Fe,O, formed on a heated steel surfac 
epitaxially on a heated macro- The tendency is for needle-like crystals to fort 
crystalline copper surface. with [110] axes normal to the surface. 
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Figure 20. General movement of carbon shot during descent of piston. 
(From Andrade and Fox (1949).) 
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Nevertheless, when bonded in a sufficiently thin layer to a steel backing, they can 
support unusually heavy loading. One such backing-bearing surface combination 
consists of a hard steel strip to which is bonded a thin copper—bronze alloy strip 
on which a layer of lead is electrodeposited. A thin ‘flash’ of indium is electro- 
deposited and later thermally diffused into the lead. The copper—bronze strip is 
necessary as, owing to the differences in crystal structure and lattice constants, it 
is difficult to deposit lead electrolytically on steel in such a way as to ensure a 
satisfactorily close co-ordination, and therefore good cohesion, between the atoms 
across the interface. ‘The total thickness of the lead—indium layer is 1/1,000 in. 
The rate of wear of both bearing surfaces is so small in relation to that of the other 
sliding surfaces of the engine that this thickness is ample (Exhibits). 

An important question relating to bearing surfaces which has so far received 
less attention than it undoubtedly merits is that of the effect of sliding on the 
surface composition of alloys. In the few cases in which this has been examined 
experimentally it has been found that polishing increases the proportion of one or 
more of the constituents in the surface of the alloy, and hence also affects the 
composition of the oxide layer. ‘Thus polishing a chromium steel alloy produces 
a surface which is either greatly enriched in chromium content or is even pure 
chromium covered by a layer of chromium oxide (Vernon, Wormwell and Nurse 
1944). <A similar effect has been observed in this laboratory with magnesium— 
aluminium alloys (De Brouckére 1945), «-brass and bronze (Bryant 1936) which 
become richer in magnesium, copper and tin respectively. Such changes may 
well be significant, not only because of their possible effect on the surface hardness 
and fusibility, but also because they must determine the nature of the oxide layer. 

The demands made on the reciprocating bearing surfaces in the cylinder of 
the internal combustion engine are in many respects much more severe than in 
the journal type of bearing. Surface temperatures are high, particularly near 
the combustion zone, and the pressure distribution on the piston crown is seldom 
uniform, especially at the start of the expansion stroke. Further, the change from 
dynamic to static conditions at the end of each stroke not only imposes unusual 
demands on the lubricant layers but also results in an uneven load distribution 
between the sliding surfaces. Finally, the reaction products cause corrosion wear 
(Williams 1937). 

Nominally, the bearing surfaces are those of the cylinder wall and piston 
rings, but the piston itself at times also rubs against the cylinder wall. For reasons 
of mechanical strength, lightness and thermal conductivity, the piston may consist 
of an aluminium alloy with copper and some magnesium, though in diesel engines 
a cast-iron skirt with an alloy steel crown is more usual. ‘The cylinder bearing 
surface is a cast-iron, austenitic steel or other ferrous alloy, sometimes with a coating 
of electrodeposited chromium. ‘The piston rings, however, are always made of 
castiron. Alloys with better mechanical properties have been tried as piston-ring 
metal in low-compression engines but their adequate lubrication presents diffi- 
culties. The graphite in the cast-iron ring helps greatly in the suppression of wear 
of both cylinder and ring (Finch, Quarrell and Wilman 1935, Finch and Whitmore 
1938). Electron diffraction has shown not only that the ‘run-in’ surfaces of the 
cylinder and cast-iron piston ring are covered with an amorphous oxide layer 
(Figure 10) but has also revealed the presence of an exceedingly thin adsorbed 
graphite layer (Figures 11 and 12), the flakes of which are orientated with their 
cleavage planes parallel to the surface. It is probable that this layer acts as a 
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temporary lubricant and protection to the oxide layer when the oil boundary 
layer breaks down.* 

The effect of direct sliding contact between the exceedingly hard corundum 
coating on the aluminium alloy piston and the oxide layer on the cylinder wall 
seems to be mitigated by the formation of magnesium spinel (Finch 1937), 
There is some evidence also that the rubbing of the corundum or spinel layers — 
against the ferric oxide layer on the cylinder wall may lead to the formation of 
iron spinels (‘Thirsk and Whitmore 1940). 


§5. THE ‘OXIDE LAYER 

In modern engineering practice, where frequent penetration of the lubricant 
boundary layer occurs, the part played by the oxide layers in facilitating sliding 
and preventing metal-to-metal contact is obviously of great importance. It is 
now known that most metals, if not indeed all, oxidize in dry air, however low the 
temperature. The primary oxide layer is formed very quickly, even at the 
temperature of liquid air, and then acts as a barrier against continued easy access 
of oxygen to the surface. Thus the subsequent rate of growth of the oxide layer 
decreases rapidly, and may become negligible at a thickness of the order of a few 
tens of atomic radii. Should further oxidation occur, it is very slow, even at room 
temperature, provided secondary effects, such as attack of the oxide films by 
adsorbed moisture, are avoided. Cabreraand Mott (1949) explain the mechanism 
of oxidation subsequent to the formation of the first oxygen layer in terms of the 
high contact potential gradient between the metal and the layer, which enables 
the metal ions to move through it without much help from thermal motion. There 
is, however, evidence which suggests that oxidation can also be due to oxygen 
gaining access to the metal through crystal boundaries. ‘Thus, we have observed 
than an anodically polished single-crystal metal surface oxidizes more slowly than 
a similarly treated microcrystalline surface. Further, on removing with citric 
acid a relatively thick oxide film grown at 180°c. on anodically polished poly- 
crystalline copper, it was seen that the crystal boundaries, barely visible before 
oxidation, had become strongly outlined and deepened (Figure 13). 

The surface oxidation of metals is, in general, a highly energetic process. 
At low temperatures, therefore, the rapid arrival of oxygen molecules and the low 
mobility of its atoms over the surface may be expected to favour the growth of 
amorphous rather than crystalline layers. There are many examples of films 
having been grown on single-crystal substrates cooled to liquid air temperatures 
and, as far as I am aware, all have been found to be amorphous or so submicro- 
scopically crystalline as to be indistinguishable from amorphous; yet the same 


* It is commonly supposed that the crystals of layer-lattice structures such as mica, graphite and 
molybdenum sulphide, which are sometimes used as dry lubricants, act like orientated packs of 
playing cards. ‘This cannot be so because, as I demonstrated earlier in this lecture, the layers can 
only be separated by cleaving in such a manner that the bonds between the layers are broken one by 
one. I also showed that, though freshly cleaved mica surfaces would adhere firmly on being first 
put together and could not be pulled apart by stress applied in a direction lying in the cleavage plane 
(i.e. attempting to rupture the bonds simultaneously), they lost this property on prolonged exposure 
to air or moisture. ‘Thus the ability of layer-lattice cleavage surfaces to slide easily over each other 
is due to adsorbed layers. _Layer-lattice crystals with their c axes parallel to, and crushed between, 
loaded sliding surfaces will be cleaved, but will also abrade the sliding surfaces. It follows that 
layer-lattice crystals in which the ratio of exposed cleavage plane area to thickness in the c direction is 
a maximum make the best form of solid lubricant. It is known from electron diffraction (Finch and 
Wilman 1936) and electron-microscopic studies (Rao, recent work in this laboratory) that the smaller 
the volume of a graphite crystal the higher is the value of the above ratio. A layer-lattice lubricant 
should, therefore, consist of crystals in the finest possible state of subdivision compatible with a 
maximum surface to thickness ratio. Sliding of the lubricated mechanism should not be relied on 
to break up the crystals, as this must result in wear. 
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substances deposited on the heated substrates have grown epitaxially. Even when 
a metal surface is oxidized at room temperature the oxide layer usually seems to be 
-amorphous. ‘This is certainly so in the cases of aluminium, chromium and iron 
(Figures 14 and 15); and recently Dr. Leu, working in my laboratory, has found 
that the oxide film formed on an anodically polished copper single-crystal face 
exposed to air at room temperature is also amorphous (Figures 16 and 17). 
On heating to above room temperature, amorphous oxide films differ widely 
‘in behaviour. Those on aluminium and chromium remain amorphous at 
_ temperatures up to about 600° c. and can grow in thickness only very slowly, if at 
all (Preston and Bircumshaw 1936, De Brouckére 1945). On copper, cuprous 
oxide grows epitaxially (Figure 18) and, at higher temperatures, becomes in turn 
superficially oxidized to cupric oxide (tenorite). These films, although coherent, 
_are fairly easily permeable by both oxygen and copper ions, and thus do not give 
a high degree of protection against oxidation. Still less protection is given by 
the crystalline iron oxide films (FeO > Fe,0,— Fe,O,) (Figure 19) formed under 
such conditions that the crystals have a marked tendency to grow outwards from 
the surface along an atom row of high atom population density. With other 
metals the crystalline oxide layer formed on heating is orientated in such a way 
that the crystals have in common a plane of high atom-population density parallel 
to the metal surface. They also cohere to each other and form a continuous 
film which hinders access of oxygen to the surface even when the metals are molten. 
Well-known examples are tin, zinc, antimony, bismuth, cadmium and lead. 
Oxide layers grown on heated metals when the rate of supply of oxygen is suffici- 
ently reduced, are often epitaxial, i.e. the crystal orientation in the oxide layer 
is related to and determined by that of the metal surface. Such an oxide layer 
grown on a single crystal is itself virtually a single crystal. Epitaxial oxide growth 
can also occur during the oxidation of polycrystalline surfaces, even though the 
metal crystals have no common direction or directions of orientation. ‘Then, 
however, the crystalline continuity of the deposited layer is further disturbed by 
the epitaxial reproduction of the irregularities which constitute the crystal 
boundaries. 

Crystalline oxide layers on the sliding surface are probably of rare occurrence. 
Although a metal surface exposed momentarily through rubbing may be extremely 
hot, it is, as we know from the formation of the Beilby layer, disturbed and very 
rapidly cooled, so that there will be a tendency for any re-formed oxide layer to be 
amorphous. Many ‘run-in’ bearing surfaces have been examined by electron 
diffraction; their oxide layers have been found to be amorphous. ‘The thickness 
of the amorphous films so far studied appears to be relatively uniform, so that the 
surface texture is a fair representation of that of the underlying metal. Indeed, 
the amorphous aluminium oxide layer has been used as a replica in the electron- 
microscopic study of the texture of aluminium surfaces. 

The adhesion of the oxide film to its substrate can be expected to be good by 
reason of its nature, though pseudomorphic strain in the initial layers may not be 
without influence (Finch and Quarrell 1933, Cabrera and Mott 1949). As the 
film is amorphous, and therefore isotropic, it should be equally strong in all 
directions inthe surface. Its cohesive strength in the case of iron and aluminium 
is shown by the way in which detached films, less than 100. thick, of amorphous 
ferric or aluminium oxide can be lifted off a water surface on a wide-meshed 
(80/inch) gauze without tearing, a process which also demonstrates their flexibility. 
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We may now examine how the oxide layer is so effective in preventing metal-to- 
metal contact when two surfaces are pressed together. When the surfaces first 
come together they make a three-point contact. This does not mean that three — 
pairs of points meet. ‘The probability is rather that three asperities of small radii _ 
of curvature will lodge on surfaces of large radii. Such a contact pair may be — 
represented by Figure 20 (Andrade and Fox 1949). If pressure is now applied, — 
the load will be distributed far more widely over the lower, flatter surface than over 
the point. Hence the point will flow and spread and, though in so doing its | 
oxide layer may be disrupted and thus expose the metal surface, the distribution 
of pressure ensured by the other, largely intact, oxide layer reduces the loading 
intensity on the opposing metal surface. The oxide layers are even more effective 
when the initial point contacts occur between areas of low curvature. On loading, 
conditions then approximate to the case of the pressing together of two plane 
surfaces separated by adhesive layers. Not until extensive plastic flow takes place — 
at both metal surfaces will their oxide films be broken and metallic contact be 
established. Since such contact cannot take place before the disruption of both 
oxide layers, it would be impossible to press two plane surfaces into metallic 
contact where one surface is soft and the other hard, unless the softer metal in _ 
flowing plastically were to shear the oxide layer off the hard surface. ‘This occurs 
when, for example, a lead column is compressed between a pair of anodically 
polished plane surfaces of polycrystalline copper; after extensive flow, the lead 
begins to form welds with the copper (Demonstration). 

We will now consider the behaviour of the oxide layers under normal sliding 
conditions, that is, after ‘running-in’ has prepared the surface for full loading. 
Where the lubricant boundary layers break down, the oxide layers must come into 
contact and slide over each other. Both films may then rupture, allowing metallic 
contact and cold or hot welding ; or one film may be ruptured and the other survive; 
or both films may survive, with or without the help of a re-forming of the lubricant 
boundary layers. Wherever they are in sliding contact, heat will be liberated in 
the oxide layers. 

‘The facts concerning the oxide layers which bear on the réle they play in 
sliding may be briefly summarized as follows. They consist of extremely thin 
amorphous films strongly held to the metal (ionic binding). The melting points 
and hardnesses of some of the oxides and other products of possible interest in 
sliding surface practice are given in the following Table. 

It will be evident from the table that the poor protection afforded by the oxide 
layers to single-crystal surfaces of copper or aluminium sliding under light loads 
at a low rate cannot be ascribed to either the refractoriness or hardness of the oxide 
layer. On the contrary, the oxide layer on the harder iron single-crystal surface 
gives a considerable measure of protection, although the fusibility and hardness of 
ferric oxide lie between those of cuprous oxide and aluminium oxide. We have — 
also seen that a thin Beilby layer on one of a pair of copper single-crystal surfaces 
rubbing together somehow increases the resistance of its oxide layer to rupture. — 

Experiments with pairs of dissimilar metals lead to similar results. Thus a 
copper single crystal sliding on aluminium under a light load almost immediately 
seizes (Demonstration), and the welds break in the body of the aluminium; 
obviously both oxide layers are injured. On the other hand, when a copper single 
crystal slides over an iron single crystal under a much heavier load, no seizure 
occurs (Demonstration). Yet the cuprous oxide film is broken in many places 


Table 1 


3 Hardness M.P. of 
Compound IMPA(ec)) (Moh) metal (° c.) 
Fe,O; (haematite) 1565 64 1535 
Fe,O, (magnetite) 1538 (decomposes) 64 
FeS, (pyrite) 1171 64 
FeS (troilite) 1193 4} 
FeCl, (laurenite) 670 Soft 
FeOCl = ae 
Fe,P 1290 = 
Fe,P 1100 == 
(Fe, Ni) Sx, (pentlandite) 878 344 
FeCr,O, (chromite) — 54 
NiO (bunsenite) 2090 54 1455 
NiS (millerite) 797 34 
(Ni, Fe)S, (bravoite) ~ 900 54-6 
CoO 1935 — 1495 
CoS (jaipurite) >1165 _- 
(Co, Ni, Fe)S. (hengleinite) ~1000 5 
er 1990 —- 1615 
CrN 1500 (decomposes) — 
Crk, 1000 — 
Ag,S (argentite) 825 24 960 
AgCl (cerargyrite) 455 14 
MnO (manganosite) 1650 54 1260 
Mn;0, (hausmannite) 1705 53-6 
MnsS, (hauerite) —- of 
Al,O; (corundum) 2050 9 660 
MoO; 700 Soft 2620 
MoS, (molybdenite) 1185 Soft 
SnO,  (cassiterite) 1127 (decomposes) 6-7 232 
SnO < 950 (decomposes) — 
SnS 882 — 
TiO, (rutile) 2130 6-64 1800 
MgO (periclase) 2500-2800 53-6 651 
ZnO (zincite) 1670 43 419 
ZnS (sphalerite) 1049 4 
Cu,0 1235 _- 1083 
CuO (tenorite) 1026 (decomposes 3-4 
to Cu,O) 

CuCl (nantokite) 422 — 
Cu.S (chalcocite) 1100 3 
CuS (covellite) 220 (decomposes) — 
PbO, (plattnerite) 290 (decomposes) 54 
PbO (litharge) 888 2 327 
PbCl, . 2PbO (mendipiie) 693 3 
PbO . PbSO, (lanarkite) 977 24 
PbSO, (anglesite) 1000 (decomposes) 3 
PbFe,O, (plumboferrite) 1530 5 
CdO 1426 3 BZ 
CdS (greenockite) 900 (sublimes) 34 
ZrO, (zirconia) >2950 63 1900 
Sb,03 656 24 630 
Sb,S, (stibnite) 550 2 
Bi,O; 820 43 Dai 
Wak 1970 — 1720 
InO, 700 (sublimes) —- 155 
Spinels (Mg, Fe, Ni, Cr, Mn, 

or ZnAl,O,) 1800 6-8 
Magnetites (Mg, Fe, Zn, 

or NiFe,O,) 1500 53-64 


Chromites (Mg or FeCr,O,) 4 


Among other surface compounds of possible interest are the selenides, tellurides, 
arsenides, phosphates and borates. 
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and the structure of the underlying metal is disturbed; the iron oxide, however, \ 
survives and thus prevents metal-to-metal contact. The oxide layers on surfacess 
of metals of high rigidity, such as hardened alloy steels, resist rupture at low sliding: 
speeds, even when heavily loaded. Examples are afforded by ball- bearing and 
gear-tooth surfaces in both of which, although the load intensity is high, the 
sliding speeds are low or very low. Such bearings, even without a lubricant, are: 
for a time capable of carrying surprisingly heavy loads at low speeds without! 
forming welds. It is evident, therefore, that the resistance of the oxide layer to) 
injury during sliding must, irrespective of its nature, depend in the first place on) 
the rigidity of the supporting metal. | 

As the severity of the sliding conditions increases, however, it becomes apparent} 
that the behaviour of the oxide layer is determined not only by the supporting: 
metal but also by its own properties. For example, when the bearing surfaces: 
consist of the same hard metal, the frequency of intermetallic contact and welding: 
increases when the sliding speed or load is raised. As we shall see later, however, , 
a change in the chemical composition of the oxide layer alone may reduce or even 
eliminate intermetallic contact and welding. 

It seems to me that these facts are best explained as follows. ‘The ratio of the: 
load to the area of contact is constant and a measure of the rigidity of the softer 
of the two metals of the opposing surfaces. Hence, with a given load, the loading; 
intensity increases with the rigidity of the softer metal or with that of the metal| 
common to both surfaces. ‘Therefore, since heat is developed in the oxide layers : 
wherever their surfaces are sliding in contact, the local temperature attained must, 
under comparable conditions, be higher with rigid metals than when one or both. 
are soft. Further, these local temperatures must increase with the speed of sliding. 
High local temperatures have, indeed, been observed during the sliding of glass 
or quartz surface pairs (Bowden, Stone and Tudor 1947), and this suggests that 
the factor limiting the possible temperature rise is the fusion of the oxide layer, 
If the softening temperature of the oxide layer is higher than that of the metal, 
the latter will melt and thus no longer support the oxide layer, and the ensuing 
hot metal-to-metal contact will lead to the formation of a hot weld. 

To test this view we may consider the hypoid gear which was introduced about 
1930. Though the sliding speeds in this type of gear are low, the load on the 
bearing surfaces is unusually high, and excessive wear due to welding was common. 
Some years previously it had been found that wear in cutting tools could be reduced 
by adding finely divided sulphur to the cutting lubricant. The sulphur seemed 
to act by forming a protective film on the tool tip. When, therefore, the best of 
the then known gear lubricants, such as castor oil and metallic soaps, had failed 
to give proper protection to the hypoid gear surfaces of hard ferrous alloy, it was a 
short but important step to the introduction by Evans (1937) of oil-soluble organic 
sulphur compounds which form surface films capable of preventing metal-to-metal 
contact under severe conditions of both loading and sliding speeds. Since then 
a bewildering variety of so-called ‘extreme-pressure’ additives has been increased 
by compounds containing chlorine or phosphorus. In 1942, Evans showed that 
the additive molecules did not perform their function by merely becoming 
adsorbed as such on the bearing surface, but found that a useful working hypothesis 
was to suppose that sulphide, chloride or phosphide layers were formed. Since 
then the formation of an unidentified iron sulphide compound, possibly FeS, 
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(Simard, Russell and Nelson 1941), of an iron phosphide (Fe,P) (Beeck, Givens 
and Williams 1941, Brummage 1942), ferrous chloride (FeCl) (Brummage 1942) 
and ferric oxychloride (FeOCI) (Brummage 1949, private communication) on 
ferrous bearing surfaces lubricated by oils containing the appropriate additives 
has been observed. Other layers produced in a similar manner have been 
identified on «-brass (cuprous chloride, CuCl) (Bryant 1936) and on silver (silver 
sulphide, Ag,S). The formation of a sulphur compound layer has also been 
confirmed by the use of an additive containing radioactive sulphur (Clark, Gallo 
and Lincoln 1943). ‘Thus Evans’ view that the extreme-pressure additive forms 
a layer by a chemical attack can be accepted. Up to the present, however, no 
reasonable hypothesis has been put forward to account for the remarkable resist- 
ance of this layer to rupture. Indeed, wear of the properly lubricated hypoid 

gear is so slight that metal-to-metal contact and welding must be of rare 
occurrence. 

Now, it is well known that an injury to the oxide layer on, for example, a ferrous 
bearing surface is almost immediately healed on renewal of the boundary layer 
owing to the action on the exposed metal of air and other oxidizing products 
dissolved in the ordinary lubricant. It is clear, therefore, that the superior 
protection afforded by the sulphide (or chloride or phosphide) layer formed by 
the extreme-pressure additive is due to some property inherent in the nature of 
this surface compound. It cannot be a question of hardness because, excepting 
pyrite, the sulphides, chlorides and phosphides of iron are softer than the oxides. 
A significant fact, however, is that, with few exceptions (the noble and platinum- 
group metals, tantalum, molybdenum and tungsten), the oxides likely to be formed 
on bearing surfaces are more refractory than the corresponding metals. On the 
other hand, the sulphides, chlorides and phosphides of the chief constituents of 
the alloy steels used in hypoid gears are more, or even much more, fusible than the 
metals. Finally, a zone of sliding contact between the oxide layers is a region of 
high temperature and, possibly, high pressure. 

With these facts in mind, it seems to me that the extreme-pressure additive 
acts in the following way. Before the frictional temperature at areas of rubbing 
contact becomes high enough to soften the underlying metal and thus rob the oxide 
layer of its support, the additive reacts with the oxide to form a layer of a com- 
pound, such as a sulphide or oxysulphide, which is more fusible than the metal. 
If the local temperature should rise sufficiently to melt the compound, some of the 
molten layer will smear and thus act for a time as a lubricant. The one or two 
molecular layers next to the metal, being more strongly held, will, however, 
remain intact and continue to resist rupture, because the melting point of the 
compound is too low to impair the rigidity of the metal and thus deprive these 
last molecular layers of their support. This view is further supported by 
Greenhill’s (1948) experiments which showed that the friction between sliding 
pairs of steel or copper surfaces was much reduced on replacing the oxide layers . 
by sulphides. 

It has already been remarked that it is a common experience that in time both 
bearing surfaces of a sliding element become polished, no matter how widely 
their melting points or hardnesses may differ. At first sight this seems to conflict 
with the results obtained by Bowden and Hughes (1937) who found that a soft 
material could polish a hard one of lower melting point without itself being polished. 
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The explanation lies in the difference between the conditions under which their 

experiments were carried out and those of normal sliding-surface practice. | 
There are, indeed, conditions in which hardness and fusibility are of secondary 
importance. For example, platinum (H.4:M.P.1,775°) is polished by a high- 

speed buff dressed with rouge (H.6}:M.P.1,565°). A diamond tool tip 

(H.10: M.P. >3,500°) is worn and smoothed by corundum (H.9:M.P.2,050°) - 
during the turning of aluminium. A rotating pad dressed with corundum 
polishes periclase (H.53-6:M.P.2,500-2,800°). Synthetic fibres wear and | 

polish the hard steel surfaces over ach ihe slide in textile machinery. Theend | 
of arod of camphor (M.P.178°c.) is caused to flow and be polished by a swiftly _ | 
moving ribbon of Wood’s metal (M.P.73°c.) (Spurr, recent experiments carried _ 
out in this laboratory). These examples suffice to show the importance of the. 
effect of a continual changing or otherwise of the rubbing surface. When two 

surfaces are in rubbing contact, the heat developed is shared by both. Ifone surface 

is always in rubbing contact, e.g. the tool tip, its temperature will exceed that of — 
the other surface which is continually changing, despite conduction losses. For 

example, a small steel specimen describing a continually changing path on a brass 

surface can be run without seizure under much more severe conditions than vice 

versa (Brummage, private communication). In both cases, the small specimen 

reaches a high temperature, which the steel, because of its higher melting point, is 

better able to withstand. This is precisely what happens on running in a sliding 

element pair consisting of a hard and a soft metal. ‘The soft metal is first polished 

and smoothed. ‘The asperities on the hard surface then form small surfaces which 

continue to rub in relatively long contact paths over rounded hill tops on the softer 

surface and are thus in their turn smoothed and polished. In general, where high 

spots are in sliding contact, the one of smaller radius of curvature will attain the 

highest temperature, and this may well exceed the melting point of the layer on 

the surface of larger radius of curvature. 


§6. THE BOUNDARY LAYER 


The function of the oxide layers is to prevent intermetallic friction with its 
more or less destructive consequences. ‘The oxide-layer (or sulphide etc.) 
component of friction, however, is still too high to be tolerated in general sliding 
practice, even supposing the layers could remain intact. It is the function of the 
boundary layers of the lubricant to replace the oxide-layer component of friction 
by one of still lower value. 

The lubricants consist, in general, of mixtures of straight-chain hydrocarbon 
molecules (paraffins) or their end-substituted products (alcohols, esters, ketones, 
acids). At temperatures below the melting point of the hydrocarbons, the chains, 
45a. apart, are grouped in bundles, like sticks in fascines. Such bundles are 
in fact crystals. When grown from the melt or solution on a solid surface the 
crystals are orientated with the chains either vertical or inclined to the surface at 
an angle characteristic of the compound. ‘The thickness of such an orientated 
crystalline layer may be as much as 50 molecules or more (for references see 
Brummage 1947), 

When such layers are heated, as they are under practical conditions of sliding, 
the fascines break up at a temperautre close to the melting point, leaving a true | 
boundary layer of molecules. The molecules attached directly to the surface 
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are normal or steeply inclined to it; above them is a layer of less pronounced 
orientation, and by the time the fourth layer is reached the orientation is practically 
non-existent (Finch and Zahoorbux 1937); that is to say, the film has now become 
fluid. ‘The maximum thickness of the boundary layer is therefore of the order of 
less than 200 a. 

The chains of the first layer in the boundary film are all adsorbed with one end 
attached to the oxide film, the radicle end in the case of end-substituted molecules. 
Owing to the polar nature of these, they are more firmly held than the normal 
parafins. Further, acids or esters may react with the oxide, sulphide or other 
layer to form soaps; in this case the relatively thick crystalline layer persists up 
to the softening temperature of the soap (Brummage 1947, Tingle 1950). 

The transition from a crystalline to a boundary film is reversible in that, 
on cooling, the crystals are re-formed (Brummage 1947). In consequence, new 
orientations of the crystals or fascines are observed after rubbing. For example, 
solidified fatty acid films are frequently found to contain crystals in which the 
molecules are aligned in the direction of sliding at an angle of about 5° to the surface 
(Finch 1938, Germer and Storks 1939). 

Boundary-layer friction is lower with end-substituted than with normal paraffins 
of similar chain lengths (Hardy 1936). This has generally been attributed to 
their stronger adhesion to the oxide layer. It is difficult, however, to see how this 
can materially reduce friction. It may be said, indeed, that present views on the 
mechanism of the facilitating of sliding by the lubricant boundary layer are 
inadequate. It is true that I have at times likened the boundary layer to a velvet 
pile or the bristles of a brush, and have also drawn attention to the thermal motion 
of the molecules of the layer and to possible effects of their polarity. But all these 
ideas are at best working hypotheses meant only to serve as guides to further 
experimental enquiry. 


§7. THE INTER-BOUNDARY LUBRICANT FILM 

The inter-boundary lubricant film exists only in regions of relatively light 
loading. Osborne Reynolds (1886) recognized its fluid character, and his theory 
of hydrodynamic flow not only represents the first broad generalization in the 
science of lubrication but has also dominated the geometrical design of the 
sliding element to this day. This is as it should be, for the shear strength of the 
fluid film is less than that between any of the other surface pairs involved in the 
sliding element. Boundary conditions of sliding should therefore be restricted 
as far as possible to those areas of heavier loading where the fluid film can no longer 
survive as such. 

In addition to keeping the boundary layers apart and reducing friction in the 
more lightly stressed regions, the fluid film has other important functions to fulfil. 
It must replenish supplies to the boundary layers, bring up oxidizing agents 
(oxygen, sulphur, halogens etc.) to repair damage sustained by the oxide layers or 
suitably to modify these, and carry off debris, clothing each particle with a protective 
boundary layer. Finally, the fluid layer also helps to carry off frictional heat. 

The main sources of wear in hydrodynamic lubrication are chemical corrosion 
and erosion. Erosive wear is chiefly due to debris rolling (Figure 6) between or 
striking sliding surfaces, and to cavitation (Munroe; see Rayleigh 1917) effects 
resulting from the collapse of bubbles suspended in the oil and in contact with 


the sliding surface. 
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§8. CONCLUSION 


The problem of the sliding element is a complex one. A distinguished 
physicist once spoke of it to me as ‘a very messy business’. This complexity 
together with the success of the empirical approach has undoubtedly deterred 
scientific research into the subject. Since Osborne Reynolds’ triumphant 
framing of the theory of hydrodynamic lubrication, only four schools have 
attempted to carry the matter further. Beilby’s (1921) work on surface flow and 
Hardy’s (1936) views on surface forces were destined, as we have seen, to prove of | 
fundamental importance. ‘The other two schools, both active to-day, are those of | 
Bowden in Cambridge and my own in South Kensington. Hitherto all these | 
workers have been concerned mainly with the metallic bearing surface and the | 
oil gua lubricant. If, in this lecture, I have laid particular emphasis on the oxide 
layer and on the oil as purveyor of oxidant to the high temperature zone of solid-to- 
solid sliding contact, it is because I believe that it is in this direction that the next 
great advance will be made. Without the oxide layer and oil-borne oxidants, 
whose interaction is confined to regions of high temperature, sliding under 
boundary conditions would be a practical impossibility. 
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ABSTRACT. It has been suggested that evaporation through a membrane or from a 
polymer film has a characteristic behaviour different from that commonly associated with a 
free liquid surface. In particular it is claimed that an increase in surface humidity produces 
an increased rate of evaporation through a membrane, though the opposite effect is well 
known for a liquid surface. This behaviour is said to be observed when diffusion within the 
membrane or film is the rate controlling process in systems for which the diffusion coefficient 
increases with increasing concentration. 

In this paper the behaviour of such diffusion-controlled evaporation processes is examined 
from the point of view of the diffusion equations and it is found that there is no theoretical 
reason to expect such processes to show anomalous behaviour. 

It is seen that when the experimental data for evaporation through a membrane are 
correctly interpreted they are in agreement with theoretical prediction. 

The so-called ‘ skin effect’ in the dry spinning of fibres is examined and appears to be 
due to complicating features other than diffusion. 


§1. INTRODUCTION 
ANY evaporative processes are diffusion-controlled in the sense that the 
rate of evaporation depends largely on the rate at which solvent or 
moisture is supplied to the evaporating surface by internal diffusion. 
The dry spinning of fibres and the loss of moisture from the body through the 
skin are examples of particular interest to the textile industry. 

It is of considerable practical importance to know how the rate of evaporation 
is influenced by external conditions, particularly those over which control can 
be exercised. ‘Thus in the evaporative (dry) spinning process, for example, 
not only the rate of production of fibres but also the shape of the filaments and 
their properties are profoundly affected by the rate of evaporation of solvent at 
different stages of the process. 

The present paper is concerned with the effect of the proportion of solvent 
vapour, or the relative humidity, in the atmosphere into which the evaporation 
takes place. In the case of evaporation from a free liquid surface it is well known 
that the higher the vapour pressure in the external atmosphere the less the rate 
of evaporation and vice versa. Correspondingly, the higher the wind speed 
past the exposed surface, the greater the rate of evaporation because the thickness 
of the boundary layer through which the vapour must pass in evaporating 
decreases with increasing wind speed. With this in mind it is somewhat sur- 
prising to find that several instances of the opposite behaviour have been reported, 
where an increase in relative humidity at the evaporating surface produces an 
increased rate of evaporation. ‘T’his is said to happen when evaporation occurs 
through an organic membrane or from a polymer film, in which cases the supply 
of solvent to the surface is governed by diffusion. ‘The following are examples 
of this anomalous behaviour. 
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Throughout the patent literature of dry spinning (see for example British 
Patent Specification 238,842) there is the suggestion that the rate of loss of 
solvent from the filaments, e.g. of acetone from cellulose acetate filaments, can 
be increased by including a large proportion of solvent vapour in the atmosphere 
through which the filaments pass after issuing from the dies. This behaviour 
is attributed to the fact that when the evaporation takes place into an atmosphere 
containing little or no solvent vapour a dry surface layer, frequently called a skin, 
is formed almost as soon as the filaments leave the dies. As it is known that the 
diffusion coefficient is very small in a near-dry polymer, the skin is considered to 
act as a barrier to the subsequent diffusion of solvent from the filament core and 
hence to slow down the evaporation in the later stages. The formation of the 
skin is prevented by using a drying atmosphere rich in solvent vapour, and it is 
therefore argued that a more rapid removal of solvent from the filaments is 
obtained. 

More recently Cassie (1949), discussing the characteristics for warmth in 
underwear fabrics, has stated that, at low temperatures, the heat loss by 
evaporation from the body increases as the relative humidity surrounding the 
skin increases. ‘This he considers to be a well-known effect for membranes such 
as the skin and in support of the observation refers to the diffusion coefficient 
for water through a keratin membrane, the coefficient being relatively very 
small at low moisture contents. 

The purpose of the present paper is to show mathematically that, if such 
anomalous behaviour exists, it cannot be attributed simply to the low value of 
the diffusion coefficient in the near-dry film or membrane. Calculations show 
that any abnormality in behaviour must be due to complicating features other 
than the simple diffusion process. In order to reduce the mathematical compli- 
cations to a minimum, the main calculations are carried out for a lamina instead 
of for a cylinder. Effects such as puckering of the surface of the filaments are 
neglected. 

Part I of the paper deals with the steady-state flow through a membrane, 
and Part II with the rate of the loss of solvent from a polymer film. 


AVR 
§2. STEADY-STATE EVAPORATION THROUGH A MEMBRANE 


(1) Theoretical 


If we have a membrane of thickness / separating a region of high relative 
humidity from one of low relative humidity, then, provided the humidities on 
the two sides are maintained constant, a steady-state transfer of moisture through 
the membrane is set up. The concentrations Cy) and C, just within the two 
surfaces of the membrane are constant and are the equilibrium regains for the 
high and low humidities respectively. A concentration gradient exists through 
the membrane, and moisture moves down the gradient by diffusion at a rate deter- 
mined by the concentrations C, and C, and the diffusion coefficient D. ‘The 
rate of evaporation F per unit area of the membrane is given by 
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Since we are concerned with a steady state, F' is the same through each section 
of the membrane, i.e. F is independent of x. Therefore we have 
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The above argument is true whether D is constant or not. In particular it is 
true when D is a function of concentration, and provided D is never negative in 
the range from C) to Cj, the integral in (4) must always increase or remain constant 
as the range of concentration is increased. ‘Thus if the high concentration Cy 
remains fixed, the rate of flow F can never decrease as a result of decreasing the 
lower concentration C,. The significance of (4) in this connection was brought 
to my notice by my colleague, Dr. G. S. Park. 
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Figure 1. Concentration distribution through a membrane for a diffusion coefficient 


increasing exponentially with concentration. Cy, C, are the concentrations on 
the high and low pressure sides respectively. 


Figure 1 shows calculated concentration distributions through the membrane 
when D is an exponential function of C such that D increases by 50-fold from 
C=0 to C=C). The concentration Cy on the high humidity side is the same 
for each curve but in one case C,=0 and in the other Cj=4C,y. The curves 
show how the concentration distribution adjusts itself so that the rate of flow 
can be greater when C)=0 even though the region of low diffusion coefficient 
is then included. Clearly a low diffusion coefficient is compensated by a high 
concentration gradient. 


(11) Experimental 
In a recent paper Cassie (1949) remarks that it is well known that an increase 
in relative humidity on the low humidity side of a membrane results in an 
increased rate of evaporation. He describes a convenient way of demonstrating 
the effect experimentally, by observing the rate of evaporation from a collodion 
bag filled with water and reports that “the evaporative loss from the bag actually 
decreases as the wind speed increases because the higher wind speed maintains 
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a low relative humidity at the surface of the bag”. The behaviour described 
appears plausible on the grounds of the intuitive reasoning outlined earlier in 
the paper, i.e. that increasing the surface humidity prevents the formation of a 
near-dry surface layer through which diffusion is slow. Thus there appears to 
be a discrepancy between the behaviour predicted by equation (4) on the one 
hand, and that expected intuitively and observed experimentally on the other. 

The original observations were made by using the collodion bag as a wet bulb 
thermometer. An increased wind speed was accompanied by an increase in 
temperature and from this it was concluded that the rate of evaporation had 
decreased. ‘The experiment has recently been repeated, however (Cassie, 
private communication), and, whilst the temperature change is confirmed, a direct 
measurement of the rate of evaporation reveals that this does in fact increase as 
the wind speed increases. This is in agreement with the mathematical prediction 
of equation (4) above. If a temperature increase accompanies the increase of 
wind speed (probably due to an increase in heat exchange from the moving air) 
it follows from (4) that the observed increase in evaporation will be greater than 
if the temperature remained constant, because the diffusion coefficient D increases 
with increasing temperature. ‘Thus we conclude that the rates of evaporation 
through a membrane and from a free liquid surface are both decreased by increasing 
the surface humidity or decreasing the wind speed. 
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§3. EVAPORATION OF SOLVENT FROM POLYMER FILMS AND 
FILAMENTS 


Physical Model on which Calculations are based 

It is assumed that the loss of solvent from the polymer film takes place by 
diffusion through the film with subsequent evaporation from the surface exposed 
to the drying air stream. On this assumption the course of drying may be calcu- 
lated from the diffusion equation. We have seen that the diffusion coefficient 
is generally concentration dependent in such systems, and the nature of this 
variation and also the law governing the surface evaporation must be known 
before the calculations can proceed. 

The calculated results shown below relate to the loss of acetone from a plane 
sheet of cellulose acetate exposed to a drying atmosphere on one surface only. 
The concentration of acetone is initially uniform throughout the sheet and constant 
temperature conditions are assumed. ‘T'wo kinds of condition at the exposed 
surface are considered: 

(i) The surface of the sheet is assumed to reach equilibrium with the drying 
atmosphere instantaneously when drying commences, i.e. if the atmosphere is 
acetone-free, the concentration at the surface falls instantaneously to zero; if 
the vapour pressure in the atmosphere here is p, the surface of the sheet 
immediately reaches the concentration which is in equilibrium with p. 

(ii) A condition expressing the rate of evaporation of acetone from the surface 
is assumed. ‘The simplest reasonable condition is that the rate of evaporation 
is directly proportional to the difference between the actual concentration, C,, 
in the surface of the sheet at any time and the concentration Cy, which would be 
in equilibrium with the vapour pressure in the atmosphere remote from the 
surface. Mathematically this means that at the surface 
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where « is the constant of proportionality and D, is the value of the diffusion 
coefficient when the concentration is C,; x is the space coordinate perpendicular 
to the face of the sheet, increasing positively inwards. 

The calculated rates of loss of solvent given below are obtained by evaluating 
solutions of the diffusion equation 


COs ce dC 
are gaia.) 5 tee at ee en ee 6 
at = (P55): 6) 


with a surface condition of the type (i) or (ii) above. Only the results of the 


calculations are given here. Details of methods of evaluating numerical solutions 
of the diffusion equation with a variable diffusion coefficient, and taking account 
of any shrinkage of the film accompanying loss of solvent, are described elsewhere 
(Crank and Henry 1949). 
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Figure 2. 


The first diffusion coefficient plotted in Figure 2 is derived from 
observations of the rate of penetration of acetone into commercial cellulose 
acetate sheet, using an interferometer technique described by Crank and 
Robinson (1950). ‘The significance of the second diffusion coefficient shown in 
Figure 2 is discussed later. 


§4. CALCULATED RATES OF DRYING AND DISTRIBUTION OF 
ACETONE 


(1) Surface taking up its final concentration instantaneously. 


This is the condition of vigorous surface evaporation. Figure 3 shows the 
calculated variation of concentration of acetone with distance through the sheet 
at three successive times for two different surface concentrations; (1) the surface 
concentration falls to zero instantaneously, (2) it falls to one-tenth of its initial 
value instantaneously. ‘There is little difference between the two curves at any 
time. In particular the total acetone content of the sheet at any time, represented 
by the area under the appropriate curve of Figure 3, is much the same whether 
the surface concentration falls to zero or to one-tenth its initial value, i.e. whether 
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a concentration range over which D is small is removed or not. In so far as there 
is a slight difference, the drier sheet at any time is obtained by putting the surface 
to zero concentration, which means vigorous evaporation into an acetone-free 
atmosphere. No advantage is gained by introducing acetone into the atmosphere. 

The reason for this is clear from the calculations, as it was in the steady-state 
case shown previously in Figure 1. When the value of D at the surface is very 
low the concentration gradient is correspondingly large so that the overall rate 
of loss of acetone, given by the product DdC/dx at the surface, does not alter 
appreciably. ‘This compensating effect is most obvious when the diffusion 
coefficient is zero over a range of low concentrations as in the second diffusion 
coeficient shown in Figure 2. This is a hypothetical coefficient chosen to 
exaggerate the effect. ‘The concentration—distance curve for this coefficient is 
sketched in Figure 3 (a) from general reasoning. The surface gradient is infinite 
over the concentration range for which D is zero (in this case for concentrations 
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Figure 3. Concentration distributions for the first diffusion coefficient of Figure 2. 


less than 0-10) so that the product can have a non-zero value. A finite gradient 
cannot develop for concentrations less than 0:10, since this would imply that 
solvent had been removed from a region of zero diffusion coefficient under the 
action of a finite gradient, which is not possible. 

Thus if we interpret the term ‘skin’ to mean a gel whose acetone content is 
less than a certain critical amount, so that the diffusion coefficient in the skin is 
low, we find that the lower the mean diffusion coefficient in the skin the greater 
the concentration gradient throughout the skin which is therefore thinner. 
With this kind of surface condition, therefore, it appears that it is not possible 
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to form, by diffusion of acetone, a skin which can seriously impede further diffusion 
of acetone. 


(11) A rate of surface evaporation proportional to C,—C). 


Figure 4 shows the concentration—distance curves calculated for the surface 
rate of loss condition (5) and the first diffusion coefficient of Figure 2. The 
curves of Figure 4 are for «J=16x10-*; in one case C)=0 and in the other 
Cy) =0:2, measured as a fraction of the initial concentration in the sheet. At 
any given time the surface is always drier for C)=0, but there is the rather 
interesting feature that the concentration—distance curves cross within the 
sheet. ‘Thus there is an interval of time during which the inside of the film is 
drier at a given time when the evaporating atmosphere contains acetone than 
when it is acetone-free. The corresponding curves for 4l=8x10-® were 
calculated, and found not to intersect. They are not sufficiently interesting 
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Figure 4. Concentration distributions for first diffusion coefficient of Figure 2, W/=16 x 10-*. 


to show. Since the curves for vigorous evaporative conditions, i.e. «= co, do 
not cross either, it is clear that there is a limited range of values of a for which 
crossing is observed. Although the magnitude of this effect will vary slightly 
with different values of «/ and Cy, and different diffusion coefficients, other 
calculations, detailed results of which are not shown, suggest that the effect will 
never be very large. Thus for the second diffusion coefficient in Figure 2 the 
separation between the concentration—distance curve is not appreciably different 
from that of the corresponding curves in Figure 4. 

It is clear from Figure 4, as from Figure 3, that the overall rate of loss of 


solvent from the film is not significantly increased by introducing acetone into 
the atmosphere. 
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§5. CYLINDRICAL FILAMENTS 
Similar calculations for the rate of loss of solvent from a cylinder show the 
same general behaviour as for the plane sheet. The most rapid overall loss of 
solvent is obtained when the drying atmosphere is solvent-free. The crossing 
of the concentration—distance curves is even less marked than in Figure 4. 


§6. CONCLUSION 

Theoretical considerations of diffusion-controlled evaporation do not appear 
to support the view that the rate of evaporation either through a membrane or 
from a polymer film can be increased by increasing the surface humidity. In 
the case of the membrane, the discrepancy between theory and experiment 
seems to have been due to a misinterpretation of the observational data and this 
discrepancy has now been removed. 

If the presence of solvent vapour in the atmosphere in other cases really does 
result in an increased rate of evaporation as reported, it appears that this must 
be due tc factors other than simple diffusion. In the case of the dry spinning 
process, for example, the tension to which the filaments are usually subjected 
whilst losing solvent may be an important complicating feature. ‘The degree of 
stretch, and hence of radial contraction, produced by a given tension must depend 
very much on whether a near-dry skin is formed or not in the early stages of 
drying. ‘Thus the main effect of solvent vapour in the atmosphere may be to allow 
the filament cross section to be decreased considerably by longitudinal stretching. 
The rate of drying is thereby increased since it is inversely proportional to the 
square of the radius. However, a detailed investigation of the cause of any ~ 
apparently anomalous behaviour is outside the scope of the present paper. In all 
probability the explanation will differ from one system to another. 
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10°, Rhodium-Platinum and Platinum, 13°% Rhodium- 
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ABSTRACT. E.M.F.—temperature reference tables are given for platinum, 10% rhodium— 
platinum and platinum, 13°% rhodium—platinum thermocouples over the range from 0° to 
1,760° c. The tables are based on the calibration of twelve thermocouples of each kind, six 
from each of two manufacturers. Each set of six thermocouples comprised two samples 
from each of ihree batches of wire. The calibrations are in terms of the International 
Temperature Scale as revised at the Ninth General Conference of Weights and Measures 
held in Paris in October 1948. 

The methods of intercomparison between the thermocouples and the various inter- 
polation instruments specified in the International Temperature Scale are described. It 
is shown that the difference between the new tables and the old are mainly due to changes in 
the International Temperature Scale, and that the new tables represent the shape of the 
E.M.F.—temperature characteristics of the thermocouples of either make sufficiently closely 
to give interpolation between calibration points with accuracy within +1°c. up to 1,063° c., 
+2° c. up to 1,550° c., and +3° c. above this temperature. 


St UNDER OD Ue EOMN 
EFERENCE tables for platinum, rhodium-platinum thermocouples of 
British origin have hitherto been prepared by the manufacturers of the 
thermocouples from skeleton tables given by the National Physical 
Laboratory in reports on the routine calibrations of sample thermocouples over a 
number of years, and from a reference table for 13°% couples published in 
1943 (Barber 1943) for the temperature range from 1,400° to 1,750°c. ‘This 
range of temperature is of importance in steel melting, and the reference table 
was prepared to meet the needs of the then recently developed quick-immersion 
technique of liquid steel. temperature measurement. Since the International 
Temperature Scale above 1,063°c. is realized by use of the optical pyrometer, 
the accurate determination of the form of the thermocouple calibration curve 
involves a comparison between thermocouple and optical pyrometer. In 1942 
other urgent work made it impossible to do this, and less direct methods had 
to be relied upon. 

A reference table of the kind here considered sets out to give the shape of 
the E.M.F.temperature curve as precisely as possible. It may be used to 
interpolate accurately between a small number of calibration points on any 
particular thermocouple, in any desired temperature range, by plotting the 
difference of the calibration points from the values in the reference table. The 
shape of the E.M.F.-temperature curve varies very slightly from one thermocouple 
to another, presumably on account of minute traces of impurity in the thermo- 
couple wires. ‘The accuracy of the reference table in representing the curvature 
of the E.M.F.temperature relation of any particular thermocouple is thus limited, 
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and the greater these variations are the more calibration points are needed to 
compose a difference curve which will give an accurate calibration of an individual 
couple. 

The tables given at the end of this paper have superseded those formerly 
issued by the makers, from which they differ quite markedly above the gold 
point (1,063°c.). This difference is mainly accounted for by the changes in 
numerical values of temperature in certain ranges resulting from the introduction, 
on Ist January 1949, of the International Temperature Scale (Comité inter- 
nationale des Poids et Mesures, 1948), revising the original 1927 scale. 


Sen DCO TO Barras) Veo GA Int ON 
The calibration of the thermocouples in terms of the International 
‘Temperature Scale involves a comparison ultimately with the interpolation 
instruments specified in the scale. For the range from 0° to 630°c. the 
platinum resistance thermometer is employed, using a quadratic equation, 


R,=R,(1 + At + BP), 


relating resistance to temperature, based on the melting point of ice (0°c.), the 

boiling point of water (100°c.), and the boiling point of sulphur (444-60° c.). 
From 630° to 1,063°c. the interpolation instrument is the platinum, 10°, 

rhodium—platinum thermocouple, employing a quadratic equation, 


e=a+bt+ct??, 


relating E.M.F. to temperature, based on the freezing points of antimony (630-5°c.), 
silver (960-8°c.) and gold (1,063°c.). 

Above 1,063°c. the optical pyrometer is used and is calibrated on the basis 
of the freezing point of gold and the Planck law of radiation with a value of 
1-438 for the constant C,. A temperature ¢ is then defined by the formula 


J, — exp{C,/A(T,+)} —1 
S55 a Exp {C,/A( Ty a tau} sa 1 ; 


where J, and /,,, are respectively the spectral energies radiated at wavelength A 
by unit area of a black body in unit time at the temperature ¢ and at the freezing 
point of gold, t,,. 7 is the melting point of ice in degrees kK. defined for the 
purpose as 273-15. 

The most significant differences in the scale outlined above compared with 
the 1927 scale are that the silver point is defined as 960-8°c. instead of 960-5°c., 
the value of C, is defined as 1-438 instead of 1-432, and the Planck Law of radiation 
is used instead of the Wien law. ‘The effect of these differences on the numerical 
values of temperature is shown in Figure 1, where it is seen that in the region 
between 630° and 1,063°c. the temperatures as defined by the 1948 scale are 
higher than those of the 1927 scale, with a maximum difference of 0:4° c. at about 
850°c., and above 1,063°c. they are lower, progressing from 0° at 1,063°c. to 
1-6° at 1,400° c. and 3-6° at 1,700°c. These changes account for the larger part 
of the departures of the E.M.F. values given in the present reference tables from 
those formerly supplied by the makers. ‘The adoption of the absolute electrical 
units in place of the international units on Ist January 1948 caused a rise in 
numerical values of voltage of 34 parts in 10°, and this change must also be taken 
into consideration when the new tables are compared with the old. 


°C. (Int 1927) minus °C. (/nt/948) 
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§3. CALIBRATION OVER THE RANGE 0° TO 1,063°c. 
The determination of the mean E.M.F.temperature curves for both the 10% 
and 13°, rhodium-platinum types of thermocouple has been based on the 


calibrations of twelve thermocouples of each kind, six being supplied by — 


Messrs. Johnson, Matthey & Co. Ltd. and six by the Shefheld Smelting Co. Ltd. 


Each set of six couples comprised two from each of three batches of wire. ‘They 


were annealed electrically at a temperature of about 1,600° c. and then calibrated 
with cold junctions in ice by direct comparison with platinum resistance thermo- 


meters over the range from 0° to 630° c. in stirred water, oil and salt baths (Grace 
and Hall 1943). 
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Figure 1. Changes in numerical value of temperature resulting from change from International 
Temperature Scale 1927 to International Temperature Scale 1948. 


The measurements at the freezing point of silver were made on an ingot of 
pure silver of about 380 gm. weight, cylindrical in shape, of 2-5 cm. diameter 
and 10cm.length. ‘The ingot was contained in a graphite crucible and the surface 
of the ingot was covered with graphite powder to a depth of about 1 cm. ‘The 
thermocouple was sheathed with a silica tube of 6 mm. external diameter and 
was immersed coaxially in the ingot to a depth of 9 cm., the silica sheath being 
protected from attack by the silver by a thin graphite sheath. Uniform 
temperature conditions over the length of the ingot were obtained by use of 
a substantial block of graphite for the crucible, which was formed by boring an 
axial hole 2-5 cm. in diameter in a cylinder about 7-7 cm. diameter and 20 cm. 
long. <A series of freezes was observed on two of the thermocouples to establish 


the validity of the point. . The reproducibility of E.M.F. of the thermocouples | 


from one freeze to another was within 1 microvolt, equivalent to better than 0-1° c. 
An upward displacement of the thermocouple in the ingot of 2 cm. during a 
freeze produced no greater change than 0-2 microvolt. Having thus established 
that satisfactory readings were being obtained, a freezing point observation was 
made on each of the 24 thermocouples in turn. 

For the gold-point observations, an ingot of 570 gm. weight, 2:5 cm. diameter 
and 6 cm. length was used. Since it is not necessary to protect the gold from the 
air as with silver, the graphite crucible was dispensed with and a silica crucible 
used with the thermocouple in a silica sheath immersed directly into the gold. 
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The crucible was placed in an axial hole bored in a nickel block in order to promote 
temperature uniformity. A similar degree of reproducibility was obtained as 
with the silver freeze. 


§4. CALIBRATIONS FROM 1,063° TO 1,760°c. 

The calibrations above the gold point involve comparison either directly 
or indirectly with the optical pyrometer. ‘The indirect method is to use the 
freezing points of pure metals already determined on the optical pyrometer 
scale. ‘The only points easily realized in the temperature range in question are 
the freezing points of palladium (1,552°c.) and platinum (1,769°c.). Other 
metals, such as nickel (freezing point 1,453°c.) and cobalt (freezing point 
1,492° c.) must be melted zm vacuo or in a reducing atmosphere. Even if all four 
of these points were used, however, a large temperature interval with no calibration 
point would be left between 1,063°c. and 1,453°c. The direct method was 
therefore adopted, and this required a black-body enclosure containing the hot 
junction of one or more thermocouples, the E.M.F. of which could be measured 
at the same time as the temperature of the enclosure was determined by means 
of the optical pyrometer. 

A diagram of the furnace used is shown in Figure 2. The black-body 
enclosure was formed by a hollow platinum cylinder 2 cm. in diameter and wall 
thickness 1 mm. The back wall of the cylinder was set at an angle to the front 
wall so that the optical pyrometer should not receive a reflected image of the 
sighting aperture, which was.a circular hole of 2 mm. diameter in the centre of 
the front wall. The furnace winding consisted of 10°, rhodium—platinum foil 
12-5 mm. wide and 0-025 mm. thick. 
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Figure 2. Central furnace tubes and black-body enclosure. 


The optical pyrometer employed was the primary standard instrument in 
which the telescope is divided to permit the use of a rotating sectored disc between 
objective and lamp. It was housed in a constant temperature room and was 
calibrated immediately before these measurements at the fundamental point 
of the scale, the gold point, using the furnace arrangement shown in Figure 3. 
The furnace tube, 2 cm. in bore, was unlagged except for a single radiation shield. 
The gold ingot, weighing about 30 gm., was contained in a steatite crucible and 
had an alumina sighting tube of 2-5 mm. bore passing through the axis. In the 
centre of this tube was a steatite plug through which was drilled a very fine hole. 
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The furnace was used in the horizontal position and, since it was symmetrical, 
could be viewed from either end. The more usual form of black-body tube 
for this purpose is closed at one end and is immersed vertically in an ingot of the 
metal. It is uncertain with this method that the pyrometer is sighted along the 
axis of the tube since there is no visible detail in the furnace when uniform 
temperature conditions exist. ‘The difficulty is overcome in the present arrange- 
ment, for the hole in the plug appears as a black spot in the field of the pyrometer 
when the alignment is correct. As there is only a small amount of furnace 
lagging, the end-cooling is small in comparison with radiation loss and, since 
the ingot is relatively short, very uniform temperature conditions are obtained. 
This is further aided by the fact that the thermal inertia of the furnace is small, 
and it is therefore only necessary to change the temperature of the windings by 
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Figure 3. Furnace and ingot assembly for gold point. 


a few degrees in order to cause the metal to melt or freeze, so that the temperature 
gradients are limited. Under these conditions well-defined melting- and freezing- 
point curves were obtained and no significant differences were detected between 
observations made by sighting into opposite ends of the furnace. 

The field of an optical pyrometer sighted into such a furnace is not uniformly 
bright. When the ingot is freezing the black-body tube appears as a bright 
central spot surrounded by darker zones and, when it is melting, as a dark spot 
with brighter surroundings. Under these conditions the matching of the 
filament against the central spot is influenced by the differences in intensity of 
the adjacent parts of the field, and in this determination led to an apparent difference 
between melts and freezes of 0-7°c. By using a supplementary lens to give a 
real image of the furnace, and at this position placing a diaphragm of slightly 
smaller diameter than the image of the black body, the field appears as a central 
bright spot in complete blackness when the pyrometer is focused on to the 
diaphragm; agreement between melts and freezes is then obtained. 

The pyrometer was then sighted into the platinum enclosure, into which 
had been inserted thermocouples calibrated at the freezing point of gold as 
already described. When the furnace was adjusted to give the gold-point 
temperature as indicated by the E.M.r. of the thermocouples, the optical pyro- 
meter readings agreed, to within 0-3° c., with those obtained by the direct method, 
thus gtving a check on the efficiency of the black-body radiator. A further 
check was obtained in a similar manner at the freezing point of palladium (1,552°c.), 
the thermocouples having been calibrated at that point by the wire method and 
the pyrometer by use of the gold point and a sectored disc of appropriate aperture. 
The mean difference between the two methods in this case was only about 0-1°c. 
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Selected thermocouples were calibrated against the optical pyrometer in 
pairs at temperatures of about 1,200°, 1,270°, 1,310°, 1,420°, 1,550° and 1,650°c., 
each of these temperatures, with the exception of the last-mentioned, being 
determined by use of a sectored disc of appropriate aperture to cut down the 
spectral brightness of the black body, at the temperature named, to that at the 
goldpoint. The reduction from 1,650° c. was achieved in two stages, 1,650°-1,200° 
and 1,200°-1,063°c. ‘The final measurement at the platinum point by the wire 
method, in which the platinum arm of the thermocouple is melted at the hot 
junction, was carried out in a manner similar to that previously described (Barber ~ 
1943). ‘The couples calibrated in this way were compared with the remainder, 
after each point, in a separate furnace held at the same temperature. At the same 
time, as a check for deterioration in E.M.F., comparison was made with a control 
couple at the gold point and no significant changes were detected. On completion 
of all the observations up to and including that at 1,650°c., the E.M.F. of the thermo- 
couples at the gold point had not changed by more than two microvolts from the 
original values. 


§5. DISCUSSION OF RESULTS 

Figure 4 shows the departures in E.M.F. of the mean of the six thermocouples 
of each make from the mean of all twelve thermocouples of one kind. ‘These 
departures nowhere exceed 5 jv. for either the 10°% or the 13°% types, and common 
reference tables may therefore be computed without significant loss in accuracy 
in representing the true shapes of the E.M.F.-temperature curves. The 
departures in E.M.F. of individual couples of the 13°{ type of each make from their 
own means are shown in Figures 5 and 6. The 10° couples showed similar 
deviations. Since the wires were not specially selected, the deviations may be 
taken as representative of those likely to be obtained with any other couples from 
the same sources. 
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Figure 4. Departures of mean of £.M.F. of six couples of each make from mean of all twelve. 


The means of the observations for the twelve 10° and for the twelve 13%, 
couples are recorded in Table 1, and from these values the reference ‘Tables 2 and 3 
have been computed by the method described in the Appendix. ‘The difference 
between the E.M.F. values given by the reference tables and the observational 
values are shown in columns 4 and 6 of Table 1. Apart from the values at 1,270° 
and 1,309°c. on the 10° couples and 1,270°c. on the 13° couples, there is no 
departure greater than 3 pv. (0:2°C.). 
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The accuracy of the measurements is estimated to be to within +0-1°c. up— 


to 1,063°c., +0:2°c. at the gold point as measured with the optical pyrometer, 
rising to +2°c. at 1,650°c. The calibration of a particular thermocouple 
obtained by measurement at a few temperatures and the use of a difference 
curve from the reference table will include these inaccuracies and any occurring 


in the calibration points. On the basis of these considerations it is estimated 


that observations at the gold and palladium points by the wire method and the 
use of the tables should yield an accuracy within +1°c. up to 1,063°c., +2°c. 
up to 1,550°c. and +3°c. above. 


‘Table 1 
Mean Observed Values of £.M.F. of twelve 13%, and twelve 10°, Thermocouples 
: 5 13°% couples 10°% couples 

Soe eae Pe aay E.M.F. Diff.* E.M.F. Ditties 
f (HV.) Cc.) (H-V.) Gren) 

Platinum wire 1769-0 21107 0-0 18692 —0:2 
( 1648-5 19514 0-0 17334 0-0 

| 1552-0 18193 0-0 16204 —— (ak 

Optical | 1420-9 16339 —0-1 14622 +0-2 
pyrometer 1308-8 14745 —0:3 13270 —0-6 
1270°3 14192 +0°6 12787 +0-6 

1201-0 13239 —(0-2 11958 0-0 

(| 1129-0 12242 +0-2 11098 +01 

Gold ingot 1063-0 11351°9 —0:1 10325°3 —0-1 
Silver ingot 960-8 9990°9 —0-2 9137-9 0-0 
630°5 5919-5 0-0 dD3955 Ori 

550-0 5009-7 0-0 4719-4 0-0 

500-0 4459-5 0-0 4219-8 0-0 

450-0 3920-7 0-0 3730°5 —0-1 

400-0 3397-1 0-0 3248-8 0-0 

350-0 2886°3 0-0 2776°6 0-0 

300-0 2391-9 0-0 2314-4 0-0 

250-0 1915-7 0-0 1865-7 0-0 

Resistance 200-0 1462-7 0-0 1434-7 0-0 
thermometer < 180-0 1288-1 +0-1 1207-7 0-0 

| 160-0 1118-7 0-0 1104-0 0:0 

140-0 954-7 0-0 944-8 0-0 

120-0 795-4 +0-1 790-9 0:0 

100-0 644-2 0-0 643-2 —0:2 

80-0 498-4 +0-1 498-9 0-0 

60-0 360-4 +0-1 362°3 +0-1 

40-0 Zits 0-0 233°3 0-0 

20:0 110-7 +0:1 1125 —0:1 

L 0-0 0-0 0-0 0-0 0-0 


* The temperature differences are (observation — reference table). 


It was mentioned earlier in this paper that the platinum, 10°% rhodium-— 
platinum thermocouple is the standard interpolation instrument for realizing the 
International ‘l'emperature Scale over the range from 630:5°c. to 1,063°c. It 
is of interest to compare the results of interpolated values obtained by the specified 
method using the quadratic equation based on the &.M.F. values at the freezing 
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points of antimony, silver and gold with those given by the interpolation methods 
described inthe Appendix. ‘Taking the values of 5,539-5, 9,137-9 and 10,325-3 pv. 
for the temperatures 630-5°, 960-8° and 1,063°c. as given in Table 1 for the 10% 
thermocouples, the following equation is derived : 


E= —315-3 +8-2305¢ +0-001674¢?, 


e mf. difference _ mécrovolt 
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Figure 5. Departures of E.M.F. of each couple from mean of six 13% couples ot make A. 
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Figure 6. Departures of E.M.F. of each couple from mean of six 13% couples of make B. 
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which gives values of E.M.F. of 6,799, 7,340 and 7,890 pv. at 750°, 800° and 850° c. 
respectively. These are in excess of the values in Table 1 by only 1, 1 and 2 pv. 
respectively. 

Another point of interest is to observe the measure of agreement given by 
the tables with the International Temperature Scale criteria for the 10% type 
of couple to establish its suitability for use as a standard instrument. ‘The 
criteria are 


Exy— Ex, =1,185 +0-158( 24, — 10,310) +3 pv. 
and 
Eyy— Egy =4,776 +0-631(E4,,— 10,310) +5 pv., 


where E,,, Ey, and Ey, are the E.M.F.’s in microvolts at the freezing points of 
gold, silver and antimony. ‘These criteria yield values on the basis of 
E,,,=10,325-3 from Table 1 of 1,187:2 and 4,784-9, and the tables give 1,185-9 
and 4,783-9, the agreement being well within the prescribed limits. 

In the introduction to this paper it was mentioned that the reference table 
for 13°% couples over the range from 1,400° to 1,750°c. was published in 1943, 
and at that time a full comparison of optical pyrometer and thermocouples could 
not be made. When this table has been corrected to the 1948 temperature scale 
and changed to absolute microvolts, the differences from the values in Table 3 
shown in Figure 7 are obtained. The 1943 table was based on thermocouples 
having a mean E.M.F. at the palladium point of 18,198 absolute microvolts. The 
corresponding value for the present tables is 18,193, so that the tables refer to 
thermocouples of approximately the same E.M.F. The differences shown in 
Figure 7 nowhere exceed the equivalent of 1° c. 
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Figure 7. Difference between present reference table and 1943 table. 
(Ref.: 3. Iron Steel Inst., No. 1, 1943, p. 205.) 
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APPENDIX 


Computation of the Reference Tables by L. Fox 
Mathematics Division, National Physical Laboratory 


There is no simple standard method for producing a smooth table, linearly 
interpolable, from a set of quantities obtained by physical measurement and 
given at unequal intervals of the argument. ‘The lack of statistical knowledge 
of the conditions of the experiments, the weights to be attached to the various 
observations, and a probable mathematical relation between the functions 
concerned suggest that graduation in the actuarial and statistical sense is here 
impossible and even superfluous. The only criterion is that the final table 
should be smooth, and that the tabular entries corresponding to the arguments 
of the given data should differ as little as possible from the observed quantities. 

The method adopted here is based ‘on this criterion, and consisted of four 
Stages: 

(1) From standard processes of interpolation, using differences and divided 
ditferences where necessary, values of the E.M.F. were obtained at equal intervals 
of 50° c. over the whole temperature range. The differences will not be smooth, 
but after a certain order will oscillate about a mean value which is approximately 
zero. Differences of orders higher than this are neglected in the interpolation 
process. In the present case the first neglected difference was the fourth. 

(ii) ‘The values so obtained at equal intervals were differenced, and again 
the differences of some order showed oscillation about zero. This difference 
was smoothed, by a kind of relaxation process, making suitable changes in the 
functional values to reduce the differences to very small numbers oscillating 
about zero. Some skill is required here to produce the desired effect with the 
minimum functional change. 

(111) The smoothed values were then interpolated to fifths on the National 
Accounting Machine, the resulting tables having an argument interval of 10° c. 
and being linearly interpolable. 

(iv) Asa final check, the temperature arguments corresponding to the observed 
voltages were obtained by linear inverse interpolation and compared with the 
experimental figures. The discrepancies are shown in Table 1, and are seen 
to lie well within the possible limits of experimental error. ‘This provides a 
complete verification of the adequacy of the process. 


The work described above was carried out in the Mathematics Division of 
the National Physical Laboratory. 
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A Joule-Thomson Cascade Liquefier for Helium 


By J. ASHMEAD 
Royal Society Mond Laboratory, Cambridge 


Communicated by the Low-Temp2rature Group, read at Group Meeting on 7th December 1949; 
MS. received 23rd January 1950. 


ABSTRACT. A helium liquefier has been constructed with an output of 3-8 litres per 
hour and a starting time of 70 minutes. The liquid hydrogen needed to cool the helium 
below its inversion temperature is provided continuously by a hydrogen liquefier within the 
same machine. The advantages of the machine are (i) short starting time, (ii) simplicity of 
operation and ease of servicing, (iii) closed hydrogen circuit. 


§1. INTRODUCTION 
EFORE the subject of the present paper is discussed, it seems desirable 
B to put on record the performance of its predecessor in the Mond Laboratory, 
the liquefier built by Kapitza (1934). Exact data for the years before the 
war are not readily available, but it is known that up to September 1939, the 
liquefier had been used on more than 350 days, and probably produced about 


2,000 litres of liquid helium. Between December 1945 and May 1949, 216 ~ 


runs were made and helium for 554 experiments was produced. ‘This is a very 
fine record for a machine which was really intended as an experimental prototype 
rather than a final model. ‘The starting time of the liquefier (using, of course, 
pre-cooling with liquid nitrogen) was only 50 minutes, and the rate of making 
liquid was 1-7 litres per hour. 

By 1947, it became clear that the larger cryostats and more ambitious research 
programmes of the post-war years made a liquefier of about twice the existing 
capacity desirable. Since careful consideration was given to the possible 
alternatives, it may be useful to record the reasons which led to the solution 
adopted, since some other laboratories may be influenced by similar considera- 
tions. ‘The new liquefier had to satisfy as many as possible of the following 
conditions: 

(i) Its output should be about 4 litres per hour. 

(11) If possible, the existing compressors and gas-holders should be used. 

(iii) The apparatus should be as safe as possible. 

(iv) The starting time should be short. 

(v) Simplicity in running and ease in servicing are desirable even at the 
expense of a slight reduction in efficiency or an increase in the consump- 
tion of liquid nitrogen. 

At first, the natural tendency was towards a re-designed liquefier of the expan- 
sion engine type. However, it appeared quite clear that no great increase in 
output was to be expected as long as the ‘through-put’ of the existing compressor 
was used, ~All that could be done was to produce a machine which might be a 
little more efficient than the Kapitza machine, and would be a good deal easier 
to service. ‘Io obtain the desired output, either a larger compressor would have 
to be installed, or else extra cooling would have to be provided by a subsidiary 
hydrogen circuit, as in the Moscow Kapitza machine. 
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In common with many other cryogenic laboratories the Mond Laboratory 
possesses a hydrogen liquefier, with a compressor and arrangements for storing 
the pure gas. It was therefore reasonable to consider at this stage those methods 
of helium liquefaction which demanded a supply of liquid hydrogen. ‘The high 
thermodynamic efficiency of a liquefier of the Simon expansion type makes it 
very attractive in those cases where all the helium required for a day’s experiment- 
ing can be produced in a single expansion, but if a number of successive com- 
pressions and expansions have to be made to provide all the liquid wanted, it 
seems less certain whether the average rate of production of liquid will be higher 
than with a continuously running, though less efficient, machine. Moreover, 
the Mond Laboratory does not possess that experience of small Simon liquefiers 
which other laboratories enjoy, so that the design of a large machine of that 
type would have been rather more of a new departure than most of the other 
possible devices. 

The other important existing type of helium liquefier is that employing 
straightforward Joule-Thomson expansion of helium pre-cooled by liquid 
hydrogen, such as that in the Kamerlingh Onnes Laboratory at Leiden. 
Examination of the performance of this machine showed that it should be easy 
to obtain the output required with the helium compressor available. On the 
other hand, such a scheme had two important disadvantages. It involved the 
use of liquid hydrogen, with its attendant dangers, and it entailed two starting 
times instead of one. ‘The second of these disadvantages could be overcome by 
making the hydrogen on a previous day and storing it (see for example Jones, 
Larsen and Simon 1948), so that the liquid helium could be produced fairly easily 
on the days when it was needed, but it seemed as though there might be another 
and more attractive possibility. 

The suggestion was that a hydrogen and helium liquefier should be combined 
within the same container. ‘The liquid hydrogen produced would never be 
decanted from the machine, but would serve only to cool the helium below its 
inversion temperature. A small liquefier of this type was described by Seiler 
(1941) but it differs in so many respects from the present machine that the two 
are hardly comparable. ‘lhe advantages of a liquefier of this type are appreciable. 
The hydrogen circuit is completely closed and the gas therefore needs purification 
only very infrequently. Further, since no liquid hydrogen ever leaves the 
apparatus, the explosion hazard is considerably reduced compared with a system 
in which liquid or gas is frequently handled. Again, the combination of the two 
liquefiers in one envelope enables a very good overall starting time to be obtained 
—certainly less than the sum of the two separate starting times. Experience with 
the existing hydrogen liquefier led us to believe that a liquefier of the Joule— 
Thomson type would be straightforward to run. ‘These arguments resulted in 
the choice of the type of liquefier described in this paper, and subsequent 
experience of the completed machine has, in our opinion, justified the decision 
made. 


§2. GENERAL DESIGN OF THE LIQUEFIER 
Before the details of the various components of the machine are examined, 
it is convenient to describe the general arrangement of the liquefier illustrated 
in the Plate. The distance from the plate near the top of the liquefier to the 
floor is about five feet. High pressure hydrogen (150 atmospheres) and high 
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pressure helium (25 atmospheres) are cooled in the top exchanger by returning: | 
hydrogen, returning helium, and nitrogen from the liquid nitrogen bottle just 


below the exchanger. The liquid nitrogen in this bottle boils at about 67°K 


under reduced pressure and the high pressure hydrogen and helium pass throuplt 
tubes soldered to the outside of the bottle and leave it at a temperature within | 
one or two degrees of its own. In the next exchanger the high pressure gases | 


are cooled by the returning gases, and after this exchanger the hydrogen expands 
at a Joule~Thomson valve to a pressure of 1-5 atmospheres. Part of the hydrogen | 
is liquefied and collects in the upper hydrogen bottle and the remainder of the | 
gas passes back up the heat interchangers. A small fraction (about one-fifth) | 
of the liquid passes through a transfer valve into the lower hydrogen bottle where > 
it boils at 6 cm. pressure, producing a temperature of 14°k. The gas from this 
lower bottle passes only through the short exchanger situated between the two 


hydrogen bottles, and then goes straight to room temperature via a wide-bore | 


tube. This causes a slight loss of thermal efficiency but solves the problem of 


producing a high pumping speed at the lower bottle. ‘Two hydrogen bottles | 


are used rather than one because if all the gas were expanded to 6 cm. pressure, 


the speed of the vacuum pump would have to be about 25 times as great. The | 


high pressure helium is cooled by passage through the tubes soldered to the 
outside of both hydrogen bottles, and by the hydrogen gas in the intermediate 
exchanger. 

The helium thus enters the final exchanger at 14°K., and after passing through 
the exchanger it expands at a Joule~Thomson valve where part of it liquefies and 
collects in the large copper bottle at the bottom of the machine. From there, the 
liquid is siphoned into the flasks through the vacuum jacketed tube which can 
be seen near the top right-hand corner of the photograph. 

An important feature of the design of this machine lies in the precautions 
taken to reduce the starting time. During the early stages of operation the 
hydrogen and helium do not return up the exchangers but pass straight from the 
bottles back to the gas-holders via ‘circulation valves’ on the top of the liquefier. 
This causes a certain waste of cold in the top exchanger, and hence a greater 
consumption of liquid nitrogen, but it results in a more rapid cooling of those 
parts of the machine below the liquid nitrogen bottle. "The wide-bore tubes to 
the circulation valves from the hydrogen and helium bottles are also connected 
via safety valves to the gas-holders in case a dangerous pressure builds up in 
either bottle. 

Some of the more detailed features of the design may now be examined. 


§3. THE HEAT INTERCHANGERS 

The top three heat exchangers call for littlecomment. They consist essentially 
of bundles of tubes soldered together. The data for their design are taken from 
Jacobs and Collins (1940), and the actual method of construction is very similar 
to that employed by Blanchard and Bittner (1942).. The bottom exchanger, on 
the other hand, is rather novel in design. At very low temperatures the high 
density of the helium allows it to flow turbulently even through very small bore 
tubes. An ideal heat exchanger would therefore consist, on the low pressure 
side, of a bundle of some fifty tubes of 1mm. bore, all in parallel. This effect can 
be produced by soldering to the high pressure tube a piece of crinkled copper 
foil, which is then wrapped round on itself and pushed inside a thin-walled 


| 
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German silver tube. Thus the spaces between the crinkles form fine bore tubes 
giving very efficient heat transfer between the gas and the copper. The copper 
is near the maximum of its thermal conductivity in this temperature range, and 
so the cold is efficiently conducted to the central tube. To prevent longitudinal 
conduction, the copper foil is cut into pieces 6 cm. long. The form eventually 
taken by this heat interchanger is rather like a multiple trombone, and the first 
few sections of it can be seen surrounding the helium bottle at the left of the 
photograph. 


§4. THE VALVES 
The Joule-Thomson valves and the transfer valve are almost identical in 
design. ‘The threaded part of the valve is at room temperature, near the top 
plate rather than close to the expansion orifice, so as to reduce the possibility of a 
valve seizing up because of impurities. The valves are so constructed that the 
whole valve, including the seating, can be withdrawn through the top plate after 
only two simple soldered joints have been broken. The design adopted for the 
circulation valves is a development from the tap used in the Kapitza liquefier 

for drawing off the liquid helium (Kapitza 1934). 


$5. THE SIPHON TUBE 
The liquid helium is drawn off from the helium bottle through a German 
silver tube over two metres long, 3 mm. in diameter and 0-1 mm. in wall thickness. 
This tube is, of course, protected by a vacuum jacket and a copper radiation shield, 
which is thermally bonded to the helium return tubes near the liquid nitrogen 
bottle. 


§6. RADIATION SHIELD 
The part of the liquefier below the nitrogen bottle is protected from the 
influx of room temperature radiation by a copper shield about 4 mm. thick, 
which is cooled by conduction from a subsidiary nitrogen bottle at its lower end. 


Sin (CANS IUINE 

The liquefier is surrounded by a brass case with a Gaco ring seal to the brass 
plate from which the liquefier hangs. A vacuum of better than 10-° mm. of 
mercury is maintained inside the casing by a diffusion pump of speed 20 litres 
per second. Even when the casing is shut off from the pump by a high vacuum 
valve, the thermal isolation is very good. After the liquefier is closed down, 
although no more liquid nitrogen is siphoned into the shield bottle, the rate of 
evaporation of helium is only 0-1 litre per hour and there is always some liquid 
hydrogen left in the hydrogen bottle the next day. 


§3. INDICATORS 

Besides the usual complement of pressure gauges, the liquefier is equipped 
with simple gas thermometers at the hydrogen and helium expansion valves. 
The levels of the liquids in the various bottles are observed by measuring on an 
inclined butyl phthalate manometer the pressure difference between the gas 
above the free surface in the bottle, and that in a tube connected to the bottom of 
the bottle. This second tube is curled in a horizontal plane just below the bottle, 
so that the free surface in the tube is always in the horizontal part of the tube, 


508 JF. Ashmead 


even if the heat flow down the tube varies appreciably. ‘The butyl phthalate | 
gauges are fitted with special self-returning safety traps, to guard against sudden | 
variations of pressure, and they are by-passed by suitable valves until conditions _ 
are steady. Despite the very small pressure differences available owing to the 


low density of liquid hydrogen, the level indicators work very satisfactorily. 


§9. GENERAL 


A great deal of attention has been paid in the design of this machine to ease of | 


servicing, and no attempt has been made to make the design compact by enclosing 


one part of the liquefier inside another. This principle has produced a rather | 


tall machine, but one in which every vital part can be easily exposed for examina- 
tion and repair. A pit has been dug in the floor under the liquefier so that the 


outer casing and the radiation shield can be removed with the minimum of | 
trouble. About 20 minutes is required to remove the case and half an hour to 
disconnect the radiation shield, and very little longer to replace them. All the | 
components hang from the top plate of the machine, which is supported by a_ 


cantilever so that everything is accessible. 


§10. PERFORMANCE 
Liquid hydrogen is produced about 30 minutes after the machine has started, 
and in a further 40 minutes the liquefier has reached equilibrium and is producing 


liquid helium. During this starting period 15 kg. of nitrogen is used in the main — 


bottle and 10 kg. to cool down the shield. Once the shield is cold, very little 
more nitrogen is needed to keep it cold. In the equilibrium state, the through-put 
of hydrogen is 21 m*/hour at a pressure of 150 atmospheres, and the helium 


through-put is 17 m?/hour at 25 atmospheres. The consumption of liquid — 
nitrogen is 8 kg/hour and the output of liquid helium is 3-8 litre/hour. If the — 
machine is operated as a hydrogen liquefier, it produces 5 litre/hour, which is — 


about 80°% of the theoretical yield. ‘The helium part of the circuit usually gives 


about 75% of the theoretical optimum, though as much as 85° has sometimes — 
been obtained. ‘The lower efficiency is accompanied by an increased pressure — 


drop, and is probably due to a slight ‘plating’ of the exchangers by residual 
impurities. 
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ABSTRACT. ‘The formal theory, based on the concept of thermal activation of domain 
processes, previously proposed to account for the phenomenon of magnetic viscosity, has 
been extended to describe time decrease of permeability. Within certain limits the pre- 
dictions of the theoretical analysis are verified by measurements made on time decrease of 
permeability using an alnico rod specimen. It is shown that the quantitative results of these 
measurements are in good agreement with those obtained from the magnetic viscosity 
experiments. ‘The possible physical mechanisms responsible for both magnetic viscosity 
and time decrease of permeability in alnico are discussed. 


§1. INTRODUCTION 


EASUREMENTS of incremental permeability p,; of a ferromagnetic 
material at any point on the magnetization curve are made by finding 
the change in magnetic induction AB produced by a known change 

in the magnetic field AH. For the majority of materials the value of p; so 
obtained is independent of the time which elapses (resting time) between the 
establishment of H and the application of the increment AH. However, when 
the incremental permeability of certain magnetic materials of low coercivity is 
examined in this way, p, is found to decrease as the resting time increases; this 
phenomenon is known as time decrease of permeability (Webb and Ford 1934), 
Snoek (1941) has shown that the effect may be observed in soft iron specimens 
only if they contain small quantities of carbon (~0-008%) or nitrogen (~0-007%) 
as impurity. He considers that the observed time decrease of permeability is 
due to a diffusion process of the impurity atoms within the iron lattice. Snoek 
has also investigated the magnetic viscosity exhibited by such soft iron specimens 
and has shown that there is an intimate relationship between the phenomena of 
magnetic viscosity and time decrease of permeability. 

In this communication a report is given of the observations which have been 
made on time decrease of permeability in alnico, a high coercivity permanent 
magnet alloy. The related phenomenon of magnetic viscosity in alnico has 
already been investigated (Street and Woolley 1949b, referred to below as [), 
and a formal theory of the phenomenon given. ‘This theory was based on the 
concept of the activation of domain processes by thermal energy. In the first 
part of the present paper, this formal theory is developed to cover the case of 
time decrease of permeability, and discussion of the underlying physical 
mechanism is deferred until a later section. 


SP AM SU DLONRNY 


As in I the underlying assumption is made that at some stage in the 
magnetization of a ferromagnetic substance a certain number of domains have 
their magnetization vectors in positions of metastable equilibrium. Hence, 
the domain vector orientation is scarcely altered by a small increase in the domain 
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energy, but a greater energy increase initiates a transition to stable equilibrium 
which, in turn, involves a discontinuous, irreversible change in vector orientation. 
Consequently, the process envisaged is essentially one of activation, and the 
critical value of the energy required to produce discontinuous rotation will be 


referred to below as the activation energy of the domain. This energy can be | 


supplied by an increase in the applied magnetic field in which case an increase 
in the magnetization of the specimen is observed. The experimental results 
given in I, which show that there is a spontaneous increase in the magnetization 
of the specimen in a steady applied field, are consistent with the assumption that 
the required activation energy may be supplied by random thermal fluctuations 
within the material. This leads to the prediction of time decrease of permeability. 
For, suppose that the applied field H (defined in §1) is of such magnitude that 
spontaneous increase of magnetization occurs, then during the resting time ¢ 
a certain number of domains are activated by thermal energy. ‘Therefore, the 
number of domains which are activated by the eventual application of AH 
decreases as a function of ¢ and, consequently, the corresponding value of AB 
decreases in a similar way, leading to an observed time decrease of permeability. 

The quantitative analysis of this process is as follows. A steady field H 
is applied to the specimen at zero time. Let N be the number of domains having 
activation energies lying between F and FE +dE at time t. 

Then, 

dN 


== — CNexp (— ir) pee (1) 


where C is a constant, k is Boltzmann’s constant and T is the absolute temperature 
ofthespecimen. In the derivation of this equation it is assumed that (a) successful 
activations do not produce additional metastable states with activation energies 
in the range from E to EF +dE and (6) each domain may be characterized by a 
constant activation energy independent of the state of the neighbouring domains, 
i.e. the magnetic interaction between domains is assumed to be negligible. 
Integrating equation (1) we have 


los == Nt-F Const, = = 20 ees (2) 


where A=Cexp(—E/kRT). 

In a general analysis, N should be taken as any arbitrary function of the 
activation energy, E. However, some definite form of this function is required 
for the further development of the analysis. In I the function was given the simple 
form defined as follows: (a) the domains have activation energies lying only in 
the range from zero to an upper limit, Zo, and (5) within this range, N at time t=0 
is independent of F and has a value pdE, where p is the distribution constant. 
As far as the phenomenon of magnetic viscosity is concerned, the result of the 
theoretical analysis using this simplification is in very good agreement with our 
experimental observations. ‘This form of the distribution is therefore retained 
in the following analysis of time decrease of permeability. Thus equation (2) 
becomes 


Nip exp (AL GES Dart ae ee (3) 


When AH is applied, the magnetization vectors of all domains having activation 
energies below some value E’ will rotate discontinuously, the actual value of E’ 
being directly determined by AH. The resultant total change in the intensity 
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of magnetization of the specimen will consist of two parts, (a) the contribution 
of the discontinuous irreversible movements, which will be denoted by AI’, 
and (b) the contribution of the reversible movements of the magnetization vectors 
of all the domains of the material. In order to determine the contribution AJ’ 
it is necessary to find the total number of domains, N,,;, having activation energies 
lying between zero and E’ at time ¢ after the field H has been established. This 
is obtained from equation (3): 


cH’ 
Nick = | P pexp(—At) dE. 


Thus, if 7 is the average value of the change in the resolved part of the domain 
magnetization along the direction of the applied field due to the activations, then 


— fe’ 
Al’ =ip / PRAY Seth a d a bor te (4) 
“0 
Changing the variable to At, equation (4) becomes 
= cA't 
Al’ = —ipkT | 5 exp — At) d(At) 
Jo At 


Be ATi Cl) — ECA 1) eee eee (5) 


where A’=Cexp(—E’/RT), and —Ei(—x) is the exponential integral, defined 
(Glaisher 1870) by 


al 
| 7 oxP (—w) du. 


The total change in the intensity of magnetization due to the application 
of AH will contain a contribution from reversible rotations of the domains, and 
this may be accounted for by the addition of a term f(AH) to the right-hand side 
of (5) which now becomes 


Ad eee PEt ACh) Fi (oo N2)| A) ee ee (6) 


In the derivation of equation (6) it has not been necessary to make any 
assumptions concerning the ranges of t and E’. A general equation of very 
similar form to that of equation (6) has been derived in the case of magnetic 
viscosity (Street and Woolley 1949 c), viz. 

AL ei phe T Eafe El = Ct) Bak Ein ok ccaacs (7) 
where A/ is the spontaneous increase in the intensity of magnetization after the 
resting time ¢, and Ay=Cexp(—£,/RT). It has been found possible to make a 
direct check of equation (7) in its most general form since it is in satisfactory 
agreement with the experimental observations made by Snoek (1947) on magnetic 


viscosity in Mn—Zn ferrite. 
A general plot of equation (7) shows that over a range of resting times ¢ 


defined by the inequalities tf) >t >t,, where 
Bina Ss Ae ] 


the value of AJ is given by zi 
Al =ipkT Int + const. 
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The convenient range of times during which experimental observations on. 
magnetic viscosity in alnico were made was found to lie well within the limits 
of t defined immediately above. ‘The experimental results were thus satisfactorily - 


represented by the preceding In? relation. 

Subject to the limitations already imposed in the derivation of equation (6), 
the latter should give a satisfactory description of time decrease of permeability, 
Since the experimental results on this phenomenon given in the present paper 
were obtained for the same alnico specimen over the same time ranges which 
were employed in the magnetic viscosity experiments, (6) is now considered 
over a similarly limited range of times, viz. ty’ >t >t,’ where 


| Ei(—Ct,’)|< | Ei(—1’t,)| | 


Se ais 9) 
ands: a Nitae te © 

Within the range of times so specified equation (6) becomes 
Alo, = —ipRT(Ind't+const.)+f(AH). ss se (10) 


The limiting values of ty’ and ¢,’ defined by (9) now depend on the value of E’, | 
and hence on the corresponding value of the applied magnetic field increment AH. © 


From the experiments on magnetic viscosity in alnico it would be expected that 
the time range from f,’ to fy’ would include the times conveniently available for 
experiment (from 30 to 3,000 sec.) when E’ is comparable with Eo, i.e. when 
relatively large values of AH are used in the experiments. 

In experimental work designed to test the predictions of equation (10) certain 
variables may be conveniently held constant. For example, if the increment 
in applied magnetic field, AH, is maintained at a constant value in a series of 
experiments, since E’ is now constant equation (10) may be written 


Al, =— pRI(nt = const.) -(AH). = eee (11) 


Again, when the time ¢ and the absolute temperature T are held constant, the 
variation of AJ, as a function of EL’ is 


AT, =ipE’ + const:4+ (AH). eee (12) 


The present experimental work was directed towards the verification of 
equations (6), (11) and (12) and the evaluation of the constants characteristic 
of the domain processes involved in the magnetization cycle. 


§3. APPARATUS AND METHOD 


The rod of alnico employed in the present investigations was that used in I 
where details of its composition, dimensions and heat treatment are given. 
‘The specimen, on which a search coil was wound, was supported symmetrically 
either in an electrically heated tubular oven, the heater of which was wound 
non-inductively and fed with D.c. or, when temperatures below that of the room 
were required, in a Dewar vessel. ‘The required external magnetic field was 


applied by means of a surrounding water-cooled solenoid capable of producing | 


a maximum field of 1,100 oersteds. ‘The current flowing in the solenoid, supplied 
by a 280 volt battery of accumulators, was measured by a potentiometer method. 
The search coil was connected to a ballistic galvanometer with a standard mutual 
inductance included in circuit for calibration purposes. ‘Two galvanometers 
were used, one of period 6 seconds and the other of period 16 seconds. 
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The temperature of the specimen was determined by a calibrated thermo- 
couple attached to its mid-point, and measurements were made at fixed temperatures 
within the range —80°c. to 130°c., the maximum temperature fluctuations about 
the steady values being approximately 2° c. 

The general principle underlying the measurements has already been outlined 
in §1. In the majority of the experiments, the steady applied field H was equal 
to the coercive field ;H, of the specimen. Initially the alnico rod was brought 
into a cyclic state of magnetization and the apparatus was then allowed to reach 
thermal equilibrium with external demagnetizing field ;H,, applied to the specimen. 
Immediately before a measurement was made one more complete hysteresis cycle 
was traversed. After the waiting period ¢ seconds had elapsed, the field increment 
AH was applied, and the resulting change in magnetic induction was measured by 
means of the ballistic galvanometer. It was found most convenient to apply the 
field increment AH to the specimen by passing a current through an auxiliary 
solenoid wound over the main solenoid. In this way it was possible to obtain 
values of AH which were consistent and reproducible to within }°%, a much 
better reproducibility than was obtained by a sudden increase of the current in the 
main magnetizing solenoid. 


84 RESULTS AND INTERPRETATION 


Measurements were first made in an attempt to verify equation (6). A 
typical set of results of variation of AJ,, with ¢ is plotted logarithmically in Figure 1. 
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AT,, is used to denote the measured value of the change of magnetization of the 
specimen. ‘Ihe measurements were made at room temperature and the increment 
AH was 6:7 oersteds for every individual determination. It was found that the 
curve drawn through these experimental points could not be represented by a 
difference of exponential integral functions as predicted by equation (6). ‘This 
apparent failure of equation (6) to describe the experimental results is explained 
by the fact that, in its derivation, only those contributions to the measured change 
of magnetization directly due to the application of AH were considered. In 
practice, however, magnetic viscosity will occur at a field strength of H+AH, 
PROC. PHYS. SOC. LXIII, 7—B 2N 
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i.e. immediately after the increment AH has been applied, and this will add a 
further contribution to the measured value of the change of induction because 
of the finite period of the ballistic galvanometer. Conclusive evidence of the 
existence of a magnetic viscosity contribution was obtained experimentally by 
making measurements on time decrease of permeability using two galvanometers 


with appreciably different periodic times, viz. 6 and 16 seconds. ‘The results of — 


measurements made at room temperature with AH =4-3 oersteds are shown in 
Figure 2. These two sets of measurements were made under identical conditions 
except for the two different galvanometers used, and it is seen that the values 
of the changes in the intensity of magnetization obtained using the 16-second 
galvanometer (curve (a@)) are greater than the corresponding values obtained 
using the 6-second period instrument. ‘Thus the longer the periodic time of 
the galvanometer the greater the contribution of the magnetic viscosity to the 
measured change in intensity of magnetization. 

In order to take account of this magnetic viscosity correction it 1s necessary 
to introduce another term into the right-hand side of equation (6). A considera- 
tion of equation (3), which gives the number of domains having activation energies 
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lying between E and E +dE as a function of E, shows the difficulty of estimating 
this correction term, ‘The form of the required distribution functions may be 
calculated from equation (3), and the results of these calculations are shown 
diagrammatically in Figure 3, where N/p is plotted as a function of E/RT for the 
three values of Ct given for the individual curves. ‘The domains represented 
by the shaded areas contribute to AJ’, whilst those represented by the unshaded 
areas are responsible for the disturbing magnetic viscosity effect. It is possible 
to obtain the analytical form of the distribution represented by the unshaded areas, 
but the resulting mathematical analysis of magnetic viscosity effects, assuming 
this form of distribution, is intractable, except for the case (a) of Figure 3, where 
the resultant distribution function at field H+AH has the simple form 
N(E)=p=constant. This corresponds to the case treated in I. 
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Because of the limitations described above, an attempt was made to eliminate 
the contribution due to magnetic viscosity from the measured value of the 
change in induction by the following experimental method. A cascade of electro- 
mechanical relays was employed to disconnect the ballistic galvanometer from 
the search coil at various short times; 5¢ (of the order of tens of milliseconds), 
after the increment AH had been applied. ‘The observed change of induction 
for constant AH and constant resting period t was found to be a continuous 
function of dt. ‘There was thus no justification for taking any particular value 
of the change of induction as corresponding to A/,,,.__ It was supposed that this 
result was due to (a) the time constant of the solenoid and battery circuit, and 
(d) the effects of eddy-current lag in the specimen, resulting in a relatively slow 
building up of the increment of field AH. Consequently, any calculation of the 
magnetic viscosity correction term which assumes an instantaneous application 
of AH is of little practical value. | 

Since the observed values of the change of induction contain an undetermined 
contribution due to magnetic viscosity, quantitative agreement between the 
experimental results and equation (6) cannot be expected. ‘This difficulty is 
avoided to a considerable extent when, in any set of measurements, the contribution 
of magnetic viscosity is constant for all the individual measurements. It is shown 
below that this condition is satisfied if, after the application of the increment AH 
at time f, the remaining domains are characterized by a constant activation energy 
distribution p. For, from equation (1), the time rate of change of the intensity 
of magnetization at t=0 due to magnetic viscosity is given by 

( 


(=) =ipC | ‘exp (—E)RT)dE 
t=0 / 0 


dt 
up xp qh . thon eel | 3) 


In I it was shown that for this alnico specimen E,)>kRT, a statement implicit 
in the first inequality of equation (8) above. ‘Therefore, to a good approximation, 
equation (13) may be written as 


dI es 
(=). = ipk LEG 


Thus the contribution to the measured change of magnetization due to 
magnetic viscosity in atime Af, determined by the time of swing of the galvanometer, 
is approximately ipkTC At, which is a constant contribution for constant p. 
Reference to Figure 3 shows that the distribution function characteristic of the 
unactivated domains remaining after the application of AH will only be constant 
if the final conditions are as suggested by the curve (a). ‘This means that for 
the maximum resting period to be used in the experiments the value of E’, and. 
consequently AH, must be of a sufficiently large magnitude. ‘Thus for these: 
cases equation (6) may be written 


AI, =ipRT[Ei( — Ct) —Ei(—2’é)] + const. +f(AH), ...... (14) 


where AJ,, denotes the value of the change of magnetization which should be 

obtained in experimental measurements. ‘The effect of magnetic viscosity is now 

included in the constant term on the right-hand side of equation (14), and now this 
2 N-2 


516 R. Street and 7. C. Woolley 


constant depends not only on the domain constants but also on the constants of the 
ballistic galvanometer and the associated circuits. The limitations of this analysis 
now become apparent, for, owing to the inclusion of the indeterminate constant 
term in equation (14), only the slopes of the experimental curves may be used 
to determine the fundamental constants of the domain processes. 

If the magnitude of AH is sufficiently large so that equation (14) holds, then, 
although the measured values of the change in magnetization should be different 
if observations are made with two ballistic galvanometers of different periodic 
times, the difference should be constant and independent of the resting period t. 
Experimental evidence for this conclusion is presented in the two sets of curves 
shown in Figures 2 and 4. In the former case, where the size of the field 
increment AH is 4:3 oersteds, the difference between the measured values with 
ballistic galvanometers of periods 6 seconds and 16 seconds is a function of the 
resting time t. However, if the experiment is repeated for a larger value of AH, 
namely 21-4 oersteds, the two sets of results differ only by a constant value. It 
will also be seen from Figure 4 that, with AH =21-4 oersteds, the inequalities (9) 
are satisfied, for a linear variation of (AJ,,, Int) is obtained. ‘The equation of this 
straight line is given by (11), after modification to take into account the contribution 
due to magnetic viscosity, i.e. 


AL = —ipkT (Int — const.) Ef AFT Const, = =e ie (15) 


From equation (15) the slope of the (AJ,,, In#) curve is —ipkT. The magni- 
tude of the slope of the lines shown in Figure 4 is 14-7 gauss, measured at a 
temperature of 14°c., which differs by less than 19%, from the value of ipkT 
derived from the previous measurements on magnetic viscosity described in I. 

The values of AJ,, were determined as a function of time, using various 
fixed values of AH. These results are plotted in Figure 5 and the values of AH 
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Figure 5. Variation of AJmas a function of the resting period t for the various values of AH (oersteds) 
indicated on the curves. Galvanometer period=6 sec. 


used are indicated on the individual curves. It is apparent from this figure that 
when AH > 18 oersteds then (a) AJ,, is proportional to Int for values of ¢ from 
30 to 600 seconds, and (4) the slopes of the individual lines are equal and have a 
common value 14:7 gauss; (temperature 14°c.). This is in agreement with 


the prediction of equation (15) that the slope of the (AJ,,, Int) curve, i.e. ipkT | 
should be independent of AH, for sufficiently large values of AH. Thus the | 
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experimentally determined slope is again in good agreement with the results of 
the magnetic viscosity measurements. ‘The family of curves in Figure 4 serves 
to define the ranges of AH and ¢ for which the inequalities given in equation (9) 
are justified. 

A further prediction of equation (15) is that, provided 7p is independent of 
temperature, the slope of the (AZ, Int) curve should be proportional to the 
absolute temperature of the specimen. Experiments were made to check this 
prediction by making measurements of AJ, as a function of resting time t when 
the specimen was maintained at temperatures of —78°, 14°, and 130°c. The 
same value of AH, 21-4 oersteds, was used in the experiments at the three different 
temperatures, and the results are plotted in Figure 6. ‘The slopes of the experi- 
mental curves have the values 9-9, 14-7 and 20-6 gauss. Thus the ratio of slope 
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Figure 6. Variation of AJm as a function of t Figure 7. Variation of Alm as a function of 
with the specimen maintained at the the change in the applied magnetic field AH 
different temperatures indicated on the for a constant resting period t=1 minute. 
curves, AH=21-4 oersteds. Galvano- Galvanometer period=6 sec. 


meter period=6 sec. 


to absolute temperature, i.e. ipk, is constant, being equal to 0-051 gauss per degree 


Centigrade in each case, to within +1°%, giving a value of ip equal to 3-7 x 1014 
gauss per erg, which is in good agreement with the value of 3-6 x 1014 gauss per erg 
derived from the measurements on magnetic viscosity given in I. 

The presence of the magnetic viscosity term precludes the possilility of 
determining the value of C from an analysis, based on equation (14), of the 
experimental results. 

The experimental verification of equation (12) is not directly possible since the 
forms of f(AH), introduced to account for reversible changes, and F(AH)=E" 
are not known ab initio. However, measurements were made of AJ, as a function 
of AA with a resting period t = 60 seconds for each measurement; the temperature 
of the specimen was constant at 20°c. The results are plotted in Figure 7 and 
it will be seen that over a range of AH from 10 to 40 oersteds, AJ,, is practically 
a linear function of AH. Qualitatively, with the constant distribution factor p 
assumed in the analysis, as indicated in Figure 3, it would be expected that AJ’, 
the contribution to the change in intensity of magnetization due to irreversible 
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activation of the domains, would reach a constant value when E’=E,. ‘This 
corresponds physically to the situation when all the irreversible domain changes 
are activated by the application of AH. ‘The consequence is that a decrease in _ 
the slope of the (AZ,,, AH) curve would be expected for values of AH above that — 
corresponding to Ey. ‘This is a possible explanation of the decrease of slope of 
the (AJ,,, AH) curve, for values of AH greater than 40 oersteds. 

It is not possible to go further in the interpretation of the results shown in 
Figure 7 in terms of equation (12), without making further assumptions. ‘The 
underlying physical mechanism responsible for the phenomena will now be 


considered. 


$5. DISCUSSION OF POSSIBLE PHYSICAL MECHANISMS 


The agreement obtained between the experimental results on time decrease 
of permeability and (a) the theoretical analysis above and (b) the experimental 
values obtained from magnetic viscosity measurements demonstrates the essential 
correctness of the activation concepts postulated to account for both phenomena. 
Any mechanism suggested must be compatible with the assumptions made in 
the development of the formal theory. One possible mechanism, already 
suggested in I, is based on the theory recently put forward by Stoner and Wohlfarth 
(1948) to account for the high values of the coercivity of permanent magnet — 
alloys. These materials are supposed to consist of small ferromagnetic islands — 
embodied in a less ferromagnetic matrix and the sizes of the islands are such as 
to preclude the formation of domain boundaries within them. It is shown that 
over certain ranges of the hysteresis. cycles of such materials magnetization 
proceeds by discontinuous and irreversible rotations of the single domain vectors. 
In a simplified case, when the islands are prolate ellipsoids orientated at random, 
the discontinuous rotations occur over a range on the hysteresis cycle lying 
between a coercivity point and saturation. At any point in this range a certain 
number of domains have their magnetization vectors in positions of metastable 
equilibrium; that is, the free energy of a domain has a minimum value both 
before and after a discontinuous change of orientation, and transition from the 
initial to the final state may only occur if sufficient energy to overcome the inter- 
vening potential barrier is supplied. ‘Thus Stoner’s analysis would account 
for the potential barriers which may be overcome by thermal activation energy. 
In addition, as has already been shown in I, this theory would explain (a) the 
assumption that during the activation process no new domains of activation 
energy & are produced by successful activation processes, and (b) the experi- 


mentally observed tact that 7p is independent of temperature. 

Néel (1949 a) has examined the phenomenon of magnetic viscosity occurring 
in materials having a small quantity of ferromagnetic substance dispersed within 
them. He supposes that the ferromagnetic substance occurs in the form of 
particles of size below that necessary for domain boundary formation, i.e. single 
domains embedded in a non-ferromagnetic matrix are present. The postulate 
of thermal activation and the resulting log? relationship obtained by Néel are 
essentially similar to those advanced by the present authors and summarized 
in a note (Street and Woolley 1949a). In addition, Néel has also suggested a 
mechanism whereby thermal deformation of the crystal lattice of the single 
domains may give rise to discontinuous rotations of the magnetization vectors 
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In a subsequent paper, Néel has considered the phenomenon of magnetic 
viscosity in ‘massive’ materials (1949 b), and has shown that displacements of 
the domain boundary walls, due to thermal activation, may occur. This is a 
possible mechanism to account for magnetic viscosity phenomena in massive 
materials, and Néel states that the general results of this analysis are similar to those 
obtained in the consideration of single domain materials. Further information 
on this work has not yet been published, but, as far as the results on alnico are 
concerned, it is difficult to see how such irreversible boundary displacements 


are consistent with a value of ip independent of temperature and the non-production 
of new metastable states after successful activations. Further, any proposed 
mechanism should be capable of giving a satisfactory account of the high coercivity 
values of alnico, and in this respect the theory of Stoner and Wohlfarth appears 
to be the more likely. 


SO shh E CT (OP HERAT TREATMENT 


Measurements of magnetic viscosity and time decrease of permeability have 
been made using alnico specimens which have been subjected to various heat 
treatments. ‘The conclusion reached is that, if the heat treatment is such as to 
produce material of highest coercivity, then the magnetic viscosity effects are 
substantially reduced. It is hoped to publish further details of this work at a 
later date. 
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ABSTRACT. Measurements have been made of the magnetic susceptibility of fairly 
pure uranium and of 5 and 10 atomic per cent alloys of iron with uranium. Using a modifi- 
cation of the Curie method, with special pole-pieces and electromagnetic dynamometer 
control, the specimens were found to be paramagnetic with no trace of ferromagnetism in the 
temperature range 20° to 350°c. At 20°c. the susceptibility of pure uranium is 
+1-:740 x 10-°x.M.U. per gm., increasing to +1:80410~-® at 350° c.; the results fit the 
equation 


Xp =32*10-U T+ 1:564.x 10 *4-24-0 x 10-/T. 


Assuming that the iron is present in form of the compound U,.Fe, it is found that the latter 
has a susceptibility of +2-036 x 10-*£.m.u. per gm. at 20° c., which remains practically 
constant with temperature. 


§1. INTRODUCTION 

N the last of a series of papers on the magnetic properties of metals dissolved 
| in mercury, referred to hereinafter as A, Bates and Somekh (1944) described 

an apparatus for measuring the magnetic susceptibility of non-homogeneous 
systems and its use in studying cadmium, gold, silver and zinc amalgams over 
wide ranges of concentrations at room temperature. In this work ferromagnetic 
impurities were rigorously avoided, so that little or no variations of field conditions 
were required. As the Curie method was used, the amalgams were enclosed in 
one of two long cylindrical tubes mounted so that they could move along magnetic 
lines of force running between special cylindrical pole pieces in such a way that 
the force exerted on each element of amalgam was as nearly as possible independent 
of its position in the tube (cf. Bates 1948). Details of the main features of the 
apparatus and the basic theory of measurement by a dynamometer method are 
given in A, pp. 183 to 185. 

It was decided, in view of its present-day importance, to measure the suscepti- 
bility of uranium, kindly supplied by the Atomic Energy Research Establishment 
at Harwell, with the above apparatus improved in certain particulars as described 
below. It was clearly necessary to deal with possible effects of ferromagnetic 
impurities and also to make measurements above room temperature. The 
specimen tube was therefore mounted inside a vertical, cylindrical electric furnace 
consisting of a copper tube wrapped with asbestos paper, wound non-inductively 
with nichrome strip, and coated with ‘purimachos’ cement. ‘The whole was 
placed inside a brass cooling jacket through which water was continuously 
circulated. A copper—constantan thermocouple with one junction placed 
directly below the specimen tube indicated the temperature; the leads to the 
junction passed through holes drilled in a ‘sindanyo’ base to which the jacket 
was fixed. 

In A, the specimen tubes were of such a length that the centre of the specimen 
was located some 18 cm. below the glass hook of the stopper, the latter being hooked 
on to the copper framework extension of the electrodynamometer with a touch of 
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paraffin wax melted upon the junction. When the furnace temperatures rose 
above 100 to 150°c. this wax joint melted and the tube sometimes changed its 
position. Accordingly, in both the glass hook and the copper hook a well defined 
V notch was filed, and a small length of plaited silk, with a loop at each end, was 
used to connect the hooks and suspend the tube. Lateral displacements due to 
temperature changes were thus effectively eliminated and the tube could hang 
vertically with the stopper in a standard position, whether empty or containing 
liquid or solid specimens which fitted the tube. Of course, solids which did 
not exactly fit the tube might displace it from the vertical and so alter the effective 
length of the couple arm of the electrodynamometer. However, by repeating 
the measurements in such a case with the tube carefully rotated through 180° so 
that the displacement is reversed, this defect could, clearly, be overcome. The 
silk connection also had the unexpected advantage that accidental collisions of 
the tube with the furnace walls did not cause trouble. 

The apparatus was calibrated with mercury distilled by the Hulett method 
in an all-glass apparatus to reduce contamination; the susceptibility of mercury 
was taken as —0-1680 x 10-* E.M.U. pergm. The specimen tubes were evacuated 
by the simple method described in A. A new electrodynamometer was constructed, 
with the coil wound on a former in order to give less distortion and to deal with 
larger couples. 

The a.c. current to the dynamometer coils was supplied by a constant voltage 
transformer delivering 230 volts r.m.s. at 500 watts, the voltage fluctuations 
being less than 0-2°%. The same electromagnet was used as in A, but it was 
supplied from a battery of heavy-duty accumulators. 


§2. EXPERIMENTAL MEASUREMENTS 

For preliminary tests a specimen of pure uranium was turned down from a 
small rod, taking the recommended health precautions during the turning and 
handling, to form a cylinder 1-5 cm. long and 5 mm. diameter weighing about 5 gm. 
To prevent rapid oxidation of the freshly turned surface the uranium was dropped 
into a dish of concentrated nitric acid—which dissolves the oxide but only attacks 
the metal slowly—dquickly transferred to a dish of absolute alcohol and washed 
thoroughly in the latter. After rapid drying with filter paper it was slipped into 
the specimen tube, which was then evacuated rapidly. Less speed in manipu- 
lation was necessary with iron—uranium alloys which oxidize much more slowly. 

The pure uranium was tested for ferromagnetism by measuring the suscepti- 
bility &,, for different fields H and plotting k, as a function of 1/H (cf. Bates 1948, 
p. 114); if no impurity is behaving as a ferromagnetic substance the graph should 
be parallel to the 1/H axis. The values of H were small enough to ensure that 
the pole-pieces were far from saturated, so that the magnetic field was always 
of the form H(dH/dr)=const. Special tests were also made to ensure that for a 
given value of H equal portions of the specimen were all exposed to the same 
magnetic force. This was done by measuring the forces on different volumes 
of liquid in the specimen tube and showing how the force depended on the volume 
occupied. Asa final result of these tests, the specimen tube was never occupied 
to a height greater than 3 cm.; the occupied portion was kept as closely as possible 
midway between the upper and lower surfaces of the pole-pieces, and a field 
current of 3-6 amp. was not exceeded. ‘This reduction in specimen length was 
compensated by using new specimen tubes of larger internal diameter. 
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For final measurements on pure uranium at room temperature a rod 1-35 cm, 


long closely fitting the specimen tube was used, and the restoring couples were | 


found for some six suitable values of the field current. For a given value of the 


field current, the two restoring couples, with the specimen tube in a chosen | 
position and in another turned through 180° with respect to the first, differed 
by less than 0:5°%. Couples on a specimen of pure mercury were measured © 


in the same way; so, too, were those on the empty specimen tube. 

Let y, and yy, be the mass susceptibilities of the specimen and Hg respectively ; 
1s lie and i) the dynamometer restoring currents for the tube containing the 
specimen of mass m,, mercury of mass my,, and vacuum, respectively, in the same 


field; then 
Xs/ Xie = (My e/Ms)[ (2? — 20°) / (Gg? — to") = Meg) LP Leta" 


A graph of 1,2 (Iy2=12, but, strictly speaking, contains a small correction for | 


twist in the suspension), plotted against J,” is shown in Figure 1, curve (a). 
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; Figure 1. Test for absence of ferromagnetism. 


It is undoubtedly linear, so that both the couple on.the uranium specimen and 
that on the mercury were strictly proportional to H(dH/dr). Consequently, 
the uranium could not have contained a ferromagnetic impurity acting as such. 
It is easy to provide a proof that ferromagnetic impurities in both the specimen 
and mercury could not give the line (a) of Figure 1. 

Applying the Gaussian theory of errors to line (a) the value of the slope 
Ty4°/12 was found, whence the final result for the susceptibility of pure uranium 
was 1-740 + 0-003 x 10° E.m.u. per gm. ‘The error +0-003 is the probable error 
to be attached to the relative measurement; the maximum error is thus +0-015 
and the standard deviation + 0-005. 

After these measurements had been completed the specimen was analysed 
and found to contain 700 parts in 10° of iron, 60 in 10° of nickel and 5 in 108 of 
cobalt. Consequently, uranium behaves like aluminium and copper in that it 
may dissolve iron in such a way that the latter cannot exhibit ferromagnetism. 
The presence of so much iron in a ferromagnetic state would have given couples 
some twenty times as large as those observed above. 
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In view of this unexpected behaviour, steps were taken to acquire specimens 
of nominal iron contents of 5 and 10 atomic per cent. The 5°%, specimen was 
too small to fit exactly into the specimen tube so it was fitted with a cylindrical 
sheath of ‘Perspex’, for which due allowance was made in the measurements. 
The results are shown in curve (6) of Figure 1; curve (c) gives the results for 
the 10%, specimen. Using the same theoretical procedure as before, the 
susceptibilities of the specimens were found to be +1-842+0-003 and 
1-924 + 0-004 x 10-® E.M.u. per gm. respectively; and, even with as much as 
10 atomic per cent impurity of iron, there was still no trace of ferromagnetism. 

The variations of the susceptibilities of the three materials with temperature, 
measured as described, are shown in Figure 2. In some of these measurements 
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Figure 2. Variation of susceptibility of uranium-iron specimens with temperature. 


convection currents caused some small displacement of the zero of the suspended 
system. Allowance for this was adequate at lower temperatures, but appeared 
uncertain above 350° c., and in view of the relatively small changes of susceptibility 
with temperature, measurements were not recorded above 350° c. 


Say DISCUSSION OF RESULTS 

It is understood that uranium forms alloys U,.Fe and U.Fe, with iron. In 
alloys of low iron content nearly all the iron would be expected to be in combination 
as U,.Fe. Assuming no interaction between uranium metal and U,.Fe the 
susceptibility of an iron—uranium alloy should be the sum of the susceptibilities 
of pure uranium and the compound U,.Fe according to Wiedemann’s law; so 
that, if a specimen contains x atomic per cent of iron, there will be present a mass 
(«My +6xM,) of U,.Fe, where My, and My, are the atomic weights of iron and 
uranium respectively. Hence, the susceptibility y, of such a preparation should 
be given by 

Xel* My, + (100 — x) My] = Xu,.1e(% Mp, + OxMy) + xy(100 — 7x) My. 
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Now the correct iron contents for the three materials as established by later 
analysis were ()-2978, 4-884 and 9-491 atomic per cent, respectively. Consequently, — 
we have the three equations: | 


23,760(1:740 +.0-003)10-8= 442-Lyy,4e+23,320x%y, vee ees (2) 
23,930(1-842 +.0-003)10-8= 7252xy, ye+15,670xy, sees (3) 
22,080(1-924 + 0-004)10- = 14,090 xu, ne + 7,991 xy. ceeeee (4) 


These three equations containing two unknowns will not give the best values 
of yy and yy, 7 if taken in pairs and treated simultaneously, as each result so 
obtained does not take account of all the experimental data. As there are only 
three equations, we used, not the method of least squares, but the more simple 
‘least sum’ method of computation. Half of equation (3) was added to each of | 
equations (2) and (4) to give two ‘normal’ equations which were solved 
simultaneously for yy and xy, giving xy = +1:740+0-003x10-° and 
Xue.Ke = +2:036 + 0-007 x 10-§ E.M.U. per gm. 

By taking values from all three curves of Figure 2 at any chosen temperature, 
then yy and xy, 7. could be obtained at that temperature by calculations exactly 
similar to those for room temperature just described. Figure 3(a\ shows the 
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Figure 3. Variation of susceptibilities of U and U,.Fe with temperature. 


curve obtained for uranium metal, and it shows that there is a steady increase 
of susceptibility with temperature of the order of 3-5 per cent over the range 
20-350°c. ‘The corresponding curve for yy,.,. is given in Figure 3(b), which 
apparently shows a very small temperature decrease of yy, », over the same range. 
The arrows at 20°c. indicate the probable error in this part of the work, viz. 
+0-007, whilst those at 200°c. show the probable error in any of the other 
temperature measurements, viz. +0-008. This error is larger, because of the 
errors to be associated with the high temperature value relative to the room 
temperature value. It is seen that the error range indicated shows that no 
great reliance can be placed on the small decrease of yy, ». with temperature, 
i.e. within probable experimental error yy,» is independent of temperature 
in the range 20-350°c. 

It must be pointed out that these results give no real evidence to support 
the assumption that the iron actually exists as U,.Fe. Now, the equation (1) 
could be modified because each specimen might be supposed to contain, for 
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example, uranium and some non-ferromagnetic form of iron not combined with 
uranium as a compound, or even to contain uranium and UFe,, which is also 
known to exist. It was therefore thought that it might be possible, by considering 
the final error to beattached to yy as derived from these various forms of equation (1), 
to obtain an indication as to which form of equation (1) is the most correct. 
Calculations were carried out using the relation 


Xs(*Mp, + (100 — x) My) = XnexMy, + xy(100 — x) My, 

i.e. assuming only uranium and iron in a paramagnetic form to be present. The 
value obtained for yy was very close to the one quoted above, and the error to 
be attached to the value was still of the order of +0-003, thus rendering the 
calculations inconclusive. 

No account has yet been taken of other ferromagnetic elements present as 
impurity, viz. cobalt and nickel. Cobalt can be completely neglected since all 
specimens contained less than 20 parts in 10°. Nickel, however, was present 


in about 95 parts in 106 and, assuming that nickel was present as U,. Ni, equation (1) 
was written in the form 


x[* My, +2 My, + (100 —« —2)My] 
=Xuere* Mp, +6xMy) + xXu..i(2Myj +62My) + xy(100 — 7% —72) My, 

where z is the atomic per cent of nickel present. From the three equations so 
obtained, one for each specimen, values of yy, xu,.7¢ and Yy,.n; could be obtained, 
and the calculations showed that, as expected, the values of yy and yy, 7. were 
not appreciably changed. ‘The error to be attached to yy, x; was, however, so 
large that the result was valueless. 

It is possible to make some deductions concerning the behaviour of pure 
uranium by noting that the experimental results of Figure 3 fit the equation 


Xy =32 x 10-UT +1-564 x 10-§ +. 24-0 x 10-8/T, 


since this shows the great preponderance of a temperature independent term in 
the expression for the susceptibility. If we isolate the last term in the above 
equation, we find that it is equivalent to a magnetic carrier of 0-21 uz obeying 
the Curie law. It is interesting that Pauling (1947, see also Zachariasen 1949), 
suggested the allocation to uranium of 5-78 instead of 6 valency electrons, from 
which we now deduce a possible magneton number of 0:22 for the electron 
residue; the agreement is somewhat striking. 

Suppose, however, that in uranium metal we postulate that the six valency 
electrons are free and that we allow —24-3 x 10-° E.M.U. per gm. atom for the 
radon core, then the gramme atomic paramagnetism at room temperature is 
438:6x10-§. Applying Stoner’s (1935, 1936) free electron theory we find that 
the theoretical susceptibility should be 32:1g/V, x 10-® where g is the number 
of free electrons per atom and JV, is the width of the energy band in electron volts. 
The experimental data give V)=0-4392 compared with the value 15-52 for com- 
pletely free electrons. Allowing for the constant term in the preceding paragraph 
the value would be approximately 46. The calculated, (Ayy/yy.AT), 
temperature decrease of susceptibility, assuming that the characteristic temperature 
of uranium is 4220°k., is —0-29 x 10-4 per degree c.; this is opposite in sign to 
the observed effect, but the discrepancy between experiment and theory may to 
some extent be attributed to the thermal expansionof the metal. Indeed, the values 
of the coefficient of linear expansion, namely 16 x 10-6 and 23 x 10-® per degree c. 
in the temperature ranges 20 to 100°c. and 300 to 400°Cc., respectively, would 
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lead us to expect a small increase of susceptibility with temperature, although 
it would certainly be much smaller than the experimental value. 

It is not thought necessary to compare our results with those of Owen (1912) 
and Honda (1910) because these workers are unlikely to have used reasonably 
pure uranium. Since this work was completed the paper by Bommer (1941) 
has come to hand. He prepared a mixture of uranium and potassium chloride 
by heating UCI, with excess metallic potassium, later distilling off the excess 
potassium. Bommer’s value of yy at room temperature is 620 x 10~® E.M.U. 
per gm. atom, and he finds it to decrease slightly with increase in temperature. 
However, as one of his specimens showed ferromagnetism at liquid air temperatures, 
while two others did not, we feel that there is no point in discussing his results in 
detail. To carry out such low temperature measurements a new apparatus is 
necessary, and this is now being made. 

Perhaps a few words are pertinent here concerning the evidence for the 
existence of 5 f electrons in uranium; much of it has recently been summarized 
by Meggers (1947). We have Mayer’s (1941) theoretical treatment of the 
potential functions of 4 f and 5 f electrons, from which he concluded that a 
series of ‘transuranic elements’ should commence filling 5 f shells at Z=90 to 92. 
We have the chemical evidence of Seaborg (1946), Emeléus (1947) and Starke 
(1947) indicating that the first 5 f electron appears at Z=90. ‘Then there is the 
spectroscopic evidence of Kiess, Humphreys and Laun (1946) who conclude from 
the spectrum of uranium that three 5 f electrons are present. The heavier 
metallic elements which show temperature increase of paramagnetism are Ru, 
Rh and Ir, belonging to the palladium and platinum groups, and it may be 
significant that in all cases the ground term is usually considered as an F state. 

To sum up, the magnitude of the susceptibility and the absence of a Curie 
or Curie-Weiss temperature relation point to the conclusion that uranium must 
be treated as a transition element. 
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ABSTRACT. In the standard form of the phase-contrast test, the light-source is a 
pinhole comparable in size with the Airy disc of the system under test. As a consequence, 
it is difficult to obtain enough light through the pinhole to carry out the test with full 
efficiency. If a slit source and a phase-changing strip are used this difficulty is overcome, 
but the properties of the test are altered; for example, the brightness distribution seen on 
the surface of a true mirror under the test is no longer radially symmetrical. 

In the present paper the properties of the phase-strip test are examined in the mathe- 
matically simplest case where the light-source is a slit of negligible width, passing mono- 
chromatic light, and it is concluded that the test can safely be used in figuring an optical 
surface to an accuracy of about one twentieth of a fringe. Below this limit, the systematic 
error inherent in the usual interpretation of the test begins to be appreciable in the case of 
slow errors, though rapid local errors are correctly seen right up to the limit of sensitivity of 
the test. 


Se iI NERODUC TITON 


HE elegant and ingenious phase-contrast method of testing the figures of 
optical surfaces, devised about sixteen years ago by F. Zernike (1934 a), 
has been described from the experimental side by C. R. Burch (1934) and 

from the theoretical by its originator and by the present writer (Zernike 1934b, 
Linfoot 1946a). Dr. Burch has described the setting up of the test essentially 
as follows: 

An artificial star of small diameter is placed at one of the stigmatic points of 
an optical system and at the other stigmatic point there is placed the small circular 
phase-retarding disc or Zernike disc. ‘The star and the Zernike disc are both 
preferably somewhat smaller than the Airy disc of the system at the points at which 
they are respectively placed. The eye of the observer is placed close behind the 
Zernike disc, and he looks either directly at the surface under test or at its image in 
any auxiliary surfaces which may come between it and him. He then sees 
interference fringes formed between the wave leaving the surface and the supple- 
mentary wave to which its passage through the Zernike disc gives rise. If the 
Zernike disc is not larger than the Airy disc produced by the system, that part of 
the supplementary wave lying within the cone of light proceeding from the system 
will be substantially spherical, so that (since the wave front leaving the system 
departs from sphericity only by the errors of the mirror under test) the fringe 
system seen in a given wavelength will be substantially the Newton-ring system 
which would be seen if a true-figure test plate were laid on, or very close to, the 
surface under test. The colour sequence of the fringes will be determined by the 
chromatic difference of phase-retardation of the Zernike disc, and can be ascer- 
tained from the colour changes which occur when the disc is moved out of focus, 
and when it is decentred. 

A practical drawback to this form of the test is the difficulty of getting enough 
light through the very small pinhole which has to be used. A high pressure 
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mercury arc used with a good-quality condenser will overcome the difficulty to 
some extent, but even so the test needs to be carried out in a darkened room, and 
in testing an unsilvered glass mirror the amount of light is inconveniently small. 

In the Foucault knife-edge test, this shortage of light can be overcome very 
simply by substituting a fine slit for the pinhole and testing with the knife-edge 
set as accurately parallel as possible to the image of the slit. The forms and 
relative intensities of the shadows seen under the test are not altered by the change, 
provided only that the slit is too short for the aberrations of the system to vary 
appreciably along its length, and the amount of light available for testing purposes 
is increased by a factor of 50 or 100. 

It was pointed out by Zernike that a slit-source can also be used in the phase- 
contrast test provided that the phase-changing disc is replaced by a phase-changing 
strip, adjusted to cover the diffraction image of the slit. Strips of suitable 
dimensions can be made without very much difficulty by pressing fine threads of 
resin between glass, or by stroking the surface of a glass plate with a ‘cutanit’ 
(molybdenum-titanium carbide) tool*, and fairly rapid irregular errors are found 
to show very nearly the same fringe-patterns under test with a strip as with a disc. 

Thus the phase-strip test shares one of the most valuable properties of the 
Zernike-disc test, that of displaying rapid irregular errors of figure on the mirror 
surface in a way which allows their size and shape to be correctly ‘seen’ by simple 
inspection of the surface under test. 

Of even greater importance to the practical mirror-maker, however, is the 
question whether the strip test can safely be used for the control of the three 
errors which are most commonly encountered during the figuring of an optical 
surface, namely primary spherical aberration, zonal error and astigmatism. ‘The 
examination of this question is the main object of the present paper. 

In §2 formulae are derived which predict the intensities seen under the test 
on a surface with arbitrary smooth errors. In §3, these formulae are applied to 
the case of a nearly true mirror and it is shown that something near to the best 
results may be expected with a quarter-wave strip, of width about one third of 
the diameter of the Airy disc. ‘The appearance, under test with such a strip, of 
small amounts of primary astigmatism, of primary spherical aberration, and of 
local zonal error is discussed in §4, and from the examination of special cases it 
is inferred that the test can safely be used for the figuring of optical surfaces to 
an accuracy of about one twentieth of a fringe. 

Throughout the analysis it is assumed that the absorption and back-reflection 
of the strip may be disregarded, that the width of the slit is small compared with 
the Airy disc of the system under test and that the light is monochromatic. In 
practice, the width of the slit would usually be comparable in size with the diameter 
of the Airy disc and the light polychromatic, but a discussion of the extent to 
which the properties of the test are modified by these circumstances is beyond 
the scope of the present paper. 


§2. DIFFRACTION THEORY OF THE PHASE-CONTRAST STRIP TEST 

We consider first the case where the light originates in a point source, the 
geometrical image of which falls on the centre of the phase-changing strip. The 
notation and assumptions are similar to those used in earlier papers (Linfoot 
1946a,b). For convenience they are given again here. 


*C.R. Burch (1934), The engraved strips advance the phase, giving ‘ positive phase contrast’. 
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All the cases wher a convergent pencil is being null-tested can be covered by 
supposing that the wave fronts originate at the surface of a nearly spherical mirror 
M which is being tested at its approximate centre of curvature. The results 
of the test are then interpreted in terms of errors of figure of the mirror M. We 
suppose that these errors of figure may amount to several wavelengths, but 
that the error slopes on the mirror are not so steep as to spread out the visible 
image to more than a moderate multiple, say 5 or 10, of the size of the Airy disc. 


Figure 1. 


‘Then it makes no appreciable difference if we suppose the wave to be leaving a 
true spherical surface M, (Figure 1), lying everywhere within a few wave- 
lengths of the surface M, and suppose that at the point (x, y, 3) on M, the 
complex displacement is 


E(x, y) = | E(x, y) | exp ((277/A)$(x, 9) 55 


thus | E(x, y)| is the amplitude at the point (x,y) on M, and ¢(x,y) the phase 
there, while A is the wavelength of the light. ‘The variation in ¢(x,y) expresses 
the distortion of the wave fronts. Outside the boundary-circle C, lying in M, 
and corresponding to the edge of the mirror, E(x, y) 1s defined as zero. 

After leaving the mirror, the light comes to a more or less imperfect focus in 
the neighbourhood of the centre of curvature O of M, at a distance s from its 
pole A; we suppose that s is large compared with the diameter of the mirror. Just 
behind O, a viewing system L receives the light and forms and image M,’ of the 
surface M; we denote by C’ the boundary of M,’. 

At O is placed the phase-changing strip which renders variations of phase in 
the wave fronts leaving M, visible as variations of intensity in the image-surface 
M,; x, and y,, measured parallel to x and y, are coordinates in the space near O. 
The central line of the strip lies along the axis Oy, in the plane T drawn through O 
at right angles to Oz. 

We define a system of coordinate-numbers (x’,y’) in the surface My’, by 
assigning to the image P’ of the point P (x,y,z) on Mg the coordinate-numbers 
x’=x,y'=y. The intensity at P’ is thus the same thing as the ‘intensity seen 
under the test’ at P. ‘To avoid unessential complications, we suppose that L 
is a perfect optical system. 

PROG. PHYS. SOC. LXIII, 7—B 20 
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By the usual arguments, based on an application of Huyghens’ principle, the _ 
complex displacement in the intermediate image surface, which we take to be the~ 
sphere S of centre A passing through O, is | 


Wu, 2) = 5 { : | 5 BOWES, 9) dd ee (2.1) 


where 


_ 2x, _ dry, 
75 Ag Fe Pon AS 


LE Sh ae (2.2) 


and the integral is only formally over an infinite domain, since the integrand | 
vanishes outside the boundary contour C. 2 

From the assumptions made concerning the form of the errors of the surface 
under test, it follows that almost all the light crosses the intermediate image surface 
in the immediate neighbourhood of the point O. 

If the phase-changing strip were not present, the complex displacement _ 
D(x’, y’) at the point (x’, y’) in the final image-surface would be obtained by | 
integrating the product 


! et. 
Fy OXP tH u+y'v)}.W(u,v) 


over the area uw? +v?<R? of the u,v plane which falls within the aperture of the 
viewing system L. (In the case of an f/10 pencil and a viewing telescope of 1 inch 
aperture, the value of R is approximately 1,300.) ‘The effect of an idealized 
Zernike strip occupying the region —c<u'<c of the intermediate image surface 
and retarding by « radians the phase of the waves passing through it can be 
represented by changing W(u,v) into e~W(u,v) throughout this region and 
leaving it unaltered elsewhere. In saying this, we assume that the effect of 
the strip is merely to change the phase of those parts of the wave which impinge 
upon it, without influencing the neighbouring parts. This assumption is never 
strictly true, but it is a permissible approximation when the pencil under test is 
of small numerical aperture. 

In reality, the phase-strip is not an infinitesimally thin lamina but a solid body 
of length 10,000A or more, of width 10 to 100A, and of thickness comparable with A. 
Now throughout that u,v region of 5S in which an alteration in the values of 
W(u, v) can cause a visible change in the appearance of the mirror under the test, 
the distance between the surface 5 and the plane T is of the order of 10-4A in a 
practical case. For example, when a 12-inch mirror is tested at f/10 in light of 
wavelength A=2 x 10-° inch, the distance between S and T is less than 5 x 10-5 
throughout the circle uw? + v?<14-7 x 104 of radius 100 times that of the Airy disc, 
while the replacement of W(u,v) by zero outside this circle does not sensibly 
affect the appearances seen under the test.* It follows that, to the order of 
approximation here involved, the idealized phase-changing lamina which replaces 
our three-dimensional phase-strip in the analysis may as well be supposed to lie 


* This may be verified experimentally by testing a nearly true mirror of long focal ratio (about 
£/50) with an iris diaphragm placed in front of the viewing telescope. (See Linfoot (1946 b) for an 
account of the corresponding experiment on changes in the appearances seen under the Foucault 
knife-edge test.) 
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in the intermediate image surface S as in the surface T. ‘The complex displace- 
ment at the point («’, y’) of the final image surface is therefore given, to a sufficient 
approximation, by the equation 


, (> eel 1 iar “ / * , 
me. y')= cag i aad exp (—7ux' —ivy’). W(u, v) dudv 


- oe (e~ *— 1) | | exp (— tux’ —ivy’). W(u, v) du dv. 


2a Yw+u<R?, —c<su<c 


The first term of (2.3) is simply the complex displacement which results when 
the mirror is imaged through the viewing system, without the interposition of 
the Zernike strip, and the second term therefore describes the effect of the strip. 

We next replace R by oo in (2.3). This amounts to neglecting the effects. 
of the finite aperture of the viewing telescope. For mirrors with errors of the 
type under consideration, the resulting change in the predicted intensities is too 
small to be detected by the eye, except at points whose distance from the rim of 
the mirror is comparable with the resolving power of the viewing system. If, in 
order to cover the possibility of a central obstruction or of asharp step on the surface 
of M, the hypotheses are widened by allowing E(x, y) to be discontinuous along a 
finite set of arcs of differentiable curves, the above assertion remains valid provided 
we replace the words ‘‘rim of the mirror”’ by ‘‘set of points of discontinuity of 
E(x,y)’. At these points the intensities depend critically on the aperture and 
optical aberrations of the viewing system, as well as on the errors of the mirror, 
and we therefore exclude them from the discussion. At the remaining points we 
obtain, on setting R= oo in (2.3), the equation 


D(x’) = = a is exp (—zux' —ivy’). W(u, v) du dv 


+ A (eis —1) | exp (—iux’ —ivy’). W(u, v) du dv. 
2a ¥ —e<u<c, —N<vV< 0 
meoaee (2.4) 


On substituting for W(u,v) from (2.1) and applying Fourier’s inversion 
formula, (2.4) becomes 
D(x’, y’) = E(x’, y’) + P(e = 1) | pau du | er B(x, y’) dx 


Pe ed Lapeer ie » Sin cia — x’ 
=Biw,y/)++(e-*-1) [ Bfx,y') AE) 


7 
J —o x—xX 


di eee (2.5) 


These transformations are justified because E(x, y) is zero everywhere outside a 
finite region of the x,y plane, is of bounded variation in this region, and is 
differentiable at the point x=x’, y=y’. 

The intensity seen under the test, at the point (x’, y’) on the mirror-surface 
or in the halo, is then measured by the quantity 


Lee Vis | LVR at) | Pa eee Oy ans soars. (2.6) 
So far, the light-source has been assumed to be a pinhole of negligible dimen- 
sions. Inthe case of aslit of negligible width, each point of the slit can be regarded 
as a point source. If each point of the primary light source were sharply focused 
2 0-2 
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on to the plane of the slit by a wide-aperture condenser system of negligible 
aberrations, the different points of the slit could be treated as independent sources. 


In such a case, the observed intensity under the test would simply be the sum of | 
the intensity contributions from the different point sources. If we suppose that — 


the length of the slit is not great enough for off-axis aberrations of the system to 
make an appreciable contribution to the apparent errors of the surface under 
test, the intensity distributions contributed by any two points of the slit differ 
only by a constant factor. Therefore, the intensity function J(«’,y’) is still 
given by equation (2.6). 

However, in a practical case the aperture of the condenser system is often 
comparable with that of the system under test and it is not true that all points of 
the primary source are sharply focused on to the plane of the slit by the condenser. 
Each luminous particle W of the primary source contributes a wave displacement 


Ew(x, y) = Ay(x, y) E(x, y) 


on the surface of My, where E(x, y) = exp {(2772/A)4(x, y)} inside C and zero outside © 
C, as before, while A(x, y) expresses the characteristics of the wave emitted by _ 


the W particle when, after passing through the condensing system and the slit, 
it arrives at the point (x, y) on the nearly spherical mirror M. Because the width 
of the slit is negligible, A(x, y) does not vary appreciably with x, and we may write, 
to a sufficient approximation 


LEG 859) = Aa YEAR) eee (2-7) 


|£w(«, y)|? measures the intensity of the wave at this point, so that if the mirror 
is to appear uniformly lit in the absence of the phase strip 


=| Ew(x, 9) P 
Ww 
must be constant over its surface. By a suitable choice of units we can write 
this condition in the form 
Ul Age =1 on Me = eee (2.8) 
Ww 
The complex displacement contribution from the W particle in the final image 


surface is then obtained on replacing E(x,y) by Aw(y)E(x,y) in the previous 
analysis; it is 


ee , (ea? 1 —in x , , =x" 
Dus") = Awl EW 9") + 5 (#1) [ An(y")BQs, 9") SAE") ay 


=Ay(y")D(x', 9’); 


where D(x’, y’) is given by (2.5). Since the waves from the different W particles 
are mutally incoherent, the total intensity seen on M, under the test is now 


2 | Dye y')P=2| daly) PID YP 


=| D(x’, 9’), 


x—x’ 


by (2.8). ‘That is, the intensities seen on M) are still given by (2.6). The same | 
is true of the intensities in the halo, since these differ from zero only in the 


y-range where (2.8) holds. 
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From (2.5) we see that the value of D(x’, y’) at any point (x’, y’) depends only 
on the figure of the mirror M along the horizontal line y=y’ and that the halo. 
is confined to a horizontal band filling the same y-range as do the points of the 
mirror. In these two respects the phase-strip test resembles the knife edge test. 
On the other hand, since the function sin {c(x—.x’)}/(~—«x’) is equivalent to a 
function of x everywhere continuous, it can be shown that there is no bright 
‘Rayleigh ring’ at the edge of the mirror, but an apparent ordinary discontinuity 
in I(x’, y’) arising from the first term on the right of (2.5). The same is true at 
a step-discontinuity on the surface of M. 


Sa NEARLY. TRUE MIRROR; CHOICE: OF PHASE STRIP 
On a nearly true mirror of uniform reflecting power, we can write 


E(x, y)=1+1E* (x,y) («?+y?<H?) 
=0 aS Lead 9 {tM ate ear Ae (3.1) 
where E*(x, y) is real and so small that its square can be neglected in comparison 


with unity, while 2H is the diameter of the mirror surface.* (2.5) then gives, 
for x? +y?<H?, 


-Y 1 — 
D(x,y) =1+iE*(x,y) —[(1 —cos«) +ising] 2 | (1 +4E*(t,y)) = ae ay, 
ay Fx 
= Wes so eas sin c(t—x) 
=1-(1—cosa)I(¥)+sina =| B*(t,y) dt 


+i| E*(,y)—sina 1(¥)—(1 —cos«) : fer E*(t,y) eae |, 


(ae 
ee (3.2) 
where Y is written for (H?—x?)'? and 
¥_ si _ Lacks : 
L(Y)= ~ iz ae di=—[Sife(Y—w)}+Sife(Y+a)}]. 1... (3.3) 
Hence 
A RAKE 
I(x, y) =| D(x, y)P=1—4sin* 5 L(Y) —L(¥)] 
: Wipes sinc(t—x 
—2sin« | LEX, y)- = | _ EMty) we ar| 
eee (3.4) 
e < 
BiG Sli Viemak*: y) hosnn (erE*(y) ao Ni 
Gy eee K 
aioe (3.5) 
where Iy(z,y)=1-4ssin? 5 1,(Y)[1—1(¥)] Bese (3.0) 


is the intensity distribution on a true mirror. 


* Thus a mirror whose errors of figure nowhere exceed 1/20 fringe may be regarded as nearly 
true in the above sense. 
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For |x|>Y, |y| <H, that is to say in the halo, we have in place of (3.2) and 


(3.5) 


sin — 


D(x, y) = —[(1—cos«) +isina] = =[ (1 +7E*(t, y)) zt dt, 


I(x, y) = | D(«,y) P=[(1 — cos)? + sin2a] 4 


x[ ores (F pee EX(1,9) Ae) ar) | 


=4 sin? = me EY ye 


SU i(a, poo ate ee ee Oe eee ee (3.7) 


‘Thus the halo is unaffected by small errors of figure on the mirror. 

A small local error on an otherwise true mirror corresponds to the case where 
E*(x,y) vanishes everywhere outside a small x interval J of length y<H. For 
each fixed y, the integral on the right of (3.5) is then O(7) times the average value 
of E*(x,y) in the interval J and the intensity distribution is therefore represented 
to a good approximation by the equation 


I(x) = T(x; v\=— 21 yj sin an Et (Ay) ee Se (3.8) 
=1- 4sin? 5 L(Y) SLY) = 21 (¥ ) sino (a yo oe (3.9) 


Equation (3.9) shows that on a small local bump the intensity is reduced by 
an amount very nearly proportional to sin« times the height of the bump. In the 
neighbourhood of such a bump the ‘background intensity’ J)(x,y) is nearly 
constant, by (3.6) and (3.3). It follows that here, as in the Zernike disc test, the 
observed intensity change is, within the limits of visual estimation, proportional 
to the height changes over the bump, and the dark patch which (with a phase- 
retarding strip) represents the error corresponds to it in shape and intensity. 
What is even more important from the practical point of view, its boundary agrees 
exactly with the boundary of the bump, so that the mirror-maker does not find 
(as he does with the Foucault test) that his polisher has inexplicably trespassed 
outside the limits of the bump.* As a result, local figuring operations can be 
carried out with ease and certainty under the phase-strip test which would be 
impracticable under the Foucault test. 

In the case where shallow irregular error covers a large part of the mirror 
surface, the approximation (3.9) is no longer applicable, and the more accurate 
formula (3.4) must be used. But since the value of the integral 


aes sin c(t —. 

(9) = = fee y) a dE ins =) eae (3.10) 
changes only slowly as w and y vary, except near the points (0, + H), the intensity 
on and near a rapid local error is still given by an equation 

I(x, y) = —2I,(Y)sinaH*(x, y) +background intensity, ...... (3: i) 


and the above conclusions remain valid, except near the top and bottom of the 
mirror disc. 


* See Ponomarev, Bull. de l Inst. Astronomique, No. 30, Leningrad 1931, and Linfoot (1948). 
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Conclusions about the sensitiveness of the test can be drawn from (3.5) and 
(3.7). By (3.7), shallow errors do not appreciably affect the halo; it follows that 
they do not appreciably alter the total amount of light seen under test on the 
mirror disc, but merely redistribute it. That is to say 


I(x, y) dx dy = {| LQry) axay =z, eos (3212) 


x+y? <H?* 


J come 


In the case of a small local bump, the second term on the right of (3.5) describes 
quantitatively the manner in which light is displaced out of the small area covered 
by the bump and the third term describes the manner in which the displaced light 
is redistributed over a horizontal band of the mirror surface. Because it is now 
smoothly spread over a much greater area than before, the intensity changes due 
to this displaced light are small in comparison with those on the bump itself. 
In the case of a ‘cobbled surface’, a large number of irregular bumps and 
hollows are distributed at random over the whole region x2 +y? <H?. Since the 
right hand side of (3.5) remains unchanged in value when a constant is added to 
E*(x, y), we may suppose that the mean of J,(Y)E*(x, y) over the surface is zero. 
‘Then it follows from (3.5) and (3.12) that 
{| AK VAN dy 0 eee Ne (3.13) 
x+y? <H* 
Since p(x,y) varies slowly in comparison with J,(Y)E*(x,y), it follows that, 
except near the points where E*(x,y) is abnormally small, and near the points 
(0, + #), the third term on the right of (3.5) is small in comparison with the second 
term. ‘That is, the resultant effect of the redistributed light is small compared 
with the intensity variations arising from the second term on the right of (3.5). 
‘The appearance of the cobbles under the phase-strip test is therefore, to a good 
approximation, a faithful representation of their true sizes and shapes, except 
near the top and bottom of the mirror surface. 

It follows from (3.4) that the intensity changes by which the errors on a nearly 
true mirror show themselves under the test are proportional to sing. ‘The 
changes are therefore greatest when the phase-strip advances or retards the phase 
by a quarter-cycle. It cannot, however, be inferred from this that the test is 
most sensitive when « = +47, since the background intensity (x, y) also depends 
on &. 

If the test is to be sensitive in detecting slow as well as rapid errors, « and c must 
be chosen so that the intensities seen under the test do not vary too much over 
the surface of atrue mirror. In particular, /)(«, vy) must not vary too much along 
the main diameter y=0, —H<x<H. Figure 2(a) shows the intensities seen 
along the main diameter and in the halo when a true mirror is tested with « = $7 
and with selected values of Hc, the phase-strip being set at focus. In practice, 
the phase-strip would be set by trial so that the disc appears as nearly as possible 
uniformly bright. Figure 2(b) shows the intensities along the main diameter at 
the focal settings which would make them as nearly uniform as possible. ‘They 
were calculated by setting E*(«, 0) =«x? in (3.5) and giving « the values, found by 
trial, —0:065, —0-13, —0-21, —0-41, —0-70 in the respective cases Hc =0-815, 
1-2, 1-4, 1:6, 2:0, 2-5. Here cis measured in uw, v units, as before, so that Hc = 3-83 

_in the case of a strip just wide enough to cover the Airy disc completely. in the 
case He=25, the above procedure gives a rather rough approximation to the 
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true values of I(x, 0). However, it is already a rather rough approximation to 
use the focal settings which minimize intensity variations along the main diameter 
instead of over the whole disc. 


He=0-815 He = 0°815 


Disk 


(a) With strip at focus. (6) With strip adjusted to give 
greatest uniformity along 
the main diameter. 


Figure 2. Intensities seen along the main diameter of a true mirror tested with quarter-wave 
phase-retarding strips of different widths. 


From Figure 2 (6) it appears that, in the case « = +47, Hc should not exceed 
about 2-0 if a true mirror is to appear uniformly illuminated along the main dia- 
meter. ‘The width of the phase-strip is then 0-52 times the diameter of the Airy 
disc. Greater uniformity along the main diameter, with only a moderate decrease 
in intensity there, is obtained by choosing Hc=1-2. ‘The appearance of the disc 
with the phase-strip set at focus is then as shown in Figure 3. There is a very 


Figure 3._Intensities seen on true mirror under test with quarter-wave phase-retarding strip 
(Hc=1°2, «=}7) set at focus. The dotted isophotes show the intensities in the halo. 
noticeable increase in brightness near the top and bottom of the disc, but its 
appearance compares favourably in uniformity with that under the Foucault test. 

A strip with «= $7 and Hc=2-0, used at the appropriate focal setting, would 
not give appreciably worse uniformity over the disc as a whole than one with 
He =1-2 and would be about equally useful in practice. However, in the critical 
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examination of weak zonal aberration something is gained by using the narrower 
strip, together with a slit mask exposing only a horizontal band of the mirror 
surface around its main diameter. By § 2, the presence of the mask does not alter 
the intensities on the exposed portion of the mirror surface. 

Increasing |«| to a value between 47 and z would, by (3.6), decrease the values. 
of I,(x,v) and accentuate their fractional variation. Some improvement in the 
sensitivity of the test for rapid error could be obtained in this way, but at the 
cost of decreasing its reliability as a means of controlling slow errors. 

We therefore select, as a typical special case on which to base the quantitative 
estimates of sensitivity in the next Section, a strip for which « =47 and Hc=1-2; 
that is to say a quarter-wave phase-retarding strip of width approximately one third. 
of the diameter of the Airy disc. With such a strip, giving ‘negative phase 
contrast’, a local bump appears darker than the neighbouring parts of the surface: 
in accordance with (3.8). 


§4. SLOW ERRORS; LOCAL ZONAL ERROR 
A slow error whose behaviour under the test is of special practical interest 
is primary astigmatism, obtained on setting 
E(w, y) =exp {ilax?+2hay +by2)} (x2 +y<HP) 
=() (eye Ei) oo antrks (4.1) 
In testing a slightly astigmatic mirror, the phase-strip is adjusted so that the 
halo appears symmetrical * and the mirror as evenly and symmetrically illumi- 
nated as possible. It results that the strip is set centrally, but that its focal setting 
depends on the orientation of the mirror and on the amount of the astigmatism. 
The effect of moving the phase-strip through a distance 6s from O, away from 
the mirror, can be found by replacing E(x, y) on the right of (2.5) by 


exp {(x* +y*) 5. B(x, ¥) 
and D(x’,y’) on the left by exp {iy(x’2+y’?)}. D(x’, y’), where 7 =7(6s/As?). 
The resulting equation gives the new value of D(x’, y’); on substituting for 
E(x, y) from (4.1) and simplifying, it becomes 


D(x’, y’) =exp {t(ax’? + 2hx'y’ + by’?)} 


ra’ + ¥ i 
«| Oe) [exp fla tn) n't +e) +2ihey") at], 
y j 


ee 


where the upper or lower alternative is taken according as x? +?= H®. 
At the focal setting 7 = —a, (4.2) gives 
Tare, yee (ays) yo Ae Tg a (4.3) 
at other focal settings this identity does not hold. Thus there is a unique focal. 
setting at which the astigmatic mirror appears skew-symmetrical under the test. 
This setting is the one which is used in practice; it will be called the preferred 
setting. (4.2) gives, at the preferred setting 


: a’ +¥ sinct . |2 
Ie’,y’)=|DG',y’)P= |, #2 | exp Zihty! “dt fee 4 (4.4) 
TES pee iY 


The integral on the right is easily evaluated in terms of the functions Si, Ci. 


* Compare (3.7) above. If the amount of astigmatism is not small, the halo is skew-- 
symmetrical by (4.3). 
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Of the coefficients a, h, b in (4.1), only h occurs in the expression (4.4) for the 
intensity at the preferred setting. It follows that a given astigmatism shows | 
itself most strongly when the mirror is oriented with its astigmatic principal axes _ 
at angles of 45° with the direction of the phase-strip. If one of these axes is 
parallel to the direction of the phase-strip, the appearance of the mirror under 
the test is the same as that of a true mirror and the astigmatism escapes — 
detection.* 
Figure 4 shows the appearance, calculated from (4.4), of + fringe of astigmatism, 
tested with a quarter-wave phase-retarding strip of width approximately one third 
-of the diameter of the Airy disc, the mirror being oriented so that its astigmatic 
principal axes make angles of 45° with the vertical strip. In this case the 


Figure 4., One tenth of a fringe of primary astigmatism (E(x, y)=exp {27ixy/10}) under 
test with a quarter-wave phase retarding strip (Hc=1:2,«=47). The dotted 
isophotes show the intensities in the halo. 


astigmatism can easily be recognized under the test, and it may be estimated 
that one-twentieth of a fringe of astigmatism should be detectable. 

Of equal practical importance is the question whether a small amount of 
primary spherical aberration can be seen and identified without difficulty. 
Figure 5 shows the appearance of one tenth of a fringe of primary spherical 
aberration, of type 


E(x, y)=exp {877 (7? —r*)}, iat ake Re (4.5) 


tested with the same strip.| From the figure it appears that a broad dark band, 
of mean radius approximately 0-7H and not much differing from a radially 
symmetrical ring in appearance, indicates correctly the position of the ‘crest 
of the doughnut’, while the brightening of the disc towards the centre and edge 
corresponds sufficiently closely with the variations in the error function 7?—74. 
Thus one tenth of a fringe of spherical over-correction presents a striking and 
easily recognized picture under the test, and it seems safe to say that one twentieth 
of a fringe could be identified without difficulty. 


* Compare the properties of astigmatism under the Foucault test, Mon. Not. R. Astr. Soc., 
1945, 105, 193. 


+ Calculated from (2.5). My thanks are due to the Cambridge Mathematical Laboratory for 
their help in computing this figure. 
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Finally, in order to form an idea of the sensitiveness of the test in detecting 
local zonal error, we consider an ‘artificial’ zone, one twentieth of a fringe high 
all over its surface, occupying the annulus 


OOH Srec0 Tae Me tae on (4.6) 


of the mirror surface. Evidently E(x, y) has step-discontinuities on its inner and 
outer boundaries of the zone, as well as at the edge of the mirror, and the im- 
mediate neighbourhood of these points, as well as of the rim of the mirror, must 
therefore be excluded from the discussion. At the remaining points, (2.5) may 
be applied with 


E(x, 9)=1 0<r<0-6H 
=e 0-6 70-7 
SF 0-7H<r<H 
=0 PSH. tee eee (4.7) 


where, as before, 7? =x?+y?. The resulting intensity distribution is represented 
in Figure 6, which shows that the zone will appear as a clearly marked dark ring 
of about half the brightness of the nearby portions of the mirror surface. 


Figure 5. One-tenth of a fringe of spherical over- Figure 6. Flat-topped local zone 
} 0:6H<r<0-7H, 


jon (E(x, 9)=exp 4 FO (=r) 
cae eS) e re 10s one twentieth of a fringe high, under test 


under test with a quarter-wave phase- 
retarding strip (Hc=1°2, «=47). The 
numbers show the relative intensities along 
the isophotes. 


with a quarter-wave phase-retarding strip 
(Hc=1'2; «=47) set at focus. The num- 
bers show the relative intensities seen in 
their immediate neighbourhoods. 


From the results of the last two sections we may infer that when carried out 
with a slit-source of negligible width the phase-strip test can safely be used to 
figure an optical surface to an accuracy of one twentieth of a fringe. Below this 
limit, the systematic errors of the test begin to reduce its usefulness for the 
control of slow errors. Rapid local errors are shown with sufficient fidelity right 
up to the limit of sensitiveness of the test, which for such errors is of the order of 
one fiftieth of a fringe. 
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LETTERS “TO TTE DEFOR 


On the Frequency Response of PbS Transistors 


As described in a previous communication (Gebbie, Banbury and Hogarth 1950), some 
specimens of (deficit type) lead sulphide have been used successfully as crystal triodes, and 
voltage gains up to 60 have now been observed. It was of interest to compare the frequency 
response of such triodes with that of germanium transistors and to examine whether this 
response could be influenced by transverse magnetic fields as observed for n-type germanium 
(Bradner Brown 1949). Figure 1 shows response curves for one Ge and three PbS triodes, 


2:0 


Output Voltage (volts) 


0-4 


| 5 10 50 100 500 ©1000 5000 
Frequency (kc/s) 


Figure 1. Frequency response of Ge and PbS transistors. 


for 4v. (1, 2, 3) and 6v. (4) on the collector. It has been found that this frequency 
response of PbS triodes is dependent on the collector voltage. ‘Thus, on a typical specimen 
for collector voltages of 1, 3, 6 and 8 volts, the corresponding frequencies at which voltage 
gain is reduced by 50% are 420, 460, 550 and 610 kc/s. Qualitatively this is expected from 
considerations of transit time of the charge carriers which pass from the emitter to the 
collector. Higher collector voltages cannot, at present, be used on PbS, but future develop- 
ments in this direction would be expected to lead to higher frequency limits. 

A magnetic field of about +3,700 gauss was applied in the manner shown in Figure 2; 


Collector Emitter 
(+ve bias) (-ve bias) 


Positive direction 
of magnetic field 


PbS 
Specimen 


Low resistance 
base electrode 


Figure 2. Direction of magnetic field relative to transistor electrodes. 


and the corresponding values of voltage gain re-determined at various frequencies. The 
results (for 6 v. on the collector) are given in the following Table, together with similar 
measurements on Ge for purposes of comparison. It should be noted that these results are 


not directly comparable with those reported by Bradner Brown (1949), whose observations 
concerned current gain. 
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Voltage Output for Fields of + 3,700, 0, —3,700 gauss. 


Frequency PbS (p-type) : 0-05 v. input Ge (n-type): 0-1 v. input 
(ke/s.) +3700 0 —3700 +3700 0 —3700 
100 1°22 1:19 1:16 0:97 1:02 1:07 
200 0-945 0:92, 0:90 0-96 1-01 1:07 
300 0-68, 0:67, 0:65, 0-87 0:93 ° 0-97 
500 URSA 0-365 0355 0-68 0-73 0:76 
700 0-255 0-24, 0:23, 0:52 0:57 0:59 
1000 0-16; 0-16; 0:16, 0:38 0-40 0-42 
1500 0-12; 0°12; 0°12, 0:24 0:26 0:28 
2000 — — — 0-17 0-18 0:20 


According to the familiar transistor concepts, the magnetic field is expected to affect the 
average path length of the charge carriers. This length should be reduced by fields which 
tend to deflect the current flow from emitter to collector towards the surface. This would 
reduce the transit time which, together with the life-time of the carriers, limits the gain and 
its frequency response. 

Although the effect of the magnetic field is small, it can be seen from its polarity that the 
charge carriers injected by the emitter into the PbS crystal must be electrons. The mech- 
anism of this triode action is expected to be closely analogous to that of the Ge p-type 
transistor as reported by Pfann and Scaff (1949). 

Thanks are due to Professor R. W. Ditchburn for placing research facilities at our 
disposal and to Dr. C. A. Hogarth for much helpful advice. 


Department of Physics, P. C. BANBURY. 
The University, H. K. HEniscu. 
Reading. 


20th April 1950. 


BraDNER Brown, C., 1949, Phys. Rev., 76, 1736. 
GeBBIE, H. A., BANBURY, P. C., and HocartH, C. A., 1950, Prac. Phys. Soc. B, 63, 371. 
Prann, W. G., and Scarr, J. H., 1949, Phys. Rev., 76, 458 (A). 


REVIEWS OF BOOKS 


Medical X-Ray Protection up to Two Mullion Volts, NATIONAL BUREAU OF 
Sranparps. Handbook 41. Pp. 41. (Washington: Superintendent of 
Documents, U.S. Department of Commerce, March 1949.) 15. 


The present handbook supersedes a previous one of the National Bureau of Standards, 
last revised in 1936. It contains general sections, rules for specific applications in various 
kilovoltage regions and particularly details of electrical protection. ‘The protective measures 
are based upon a permissible dosage-rate of 0:3r/week in place of 0:5 r/week adopted in the 
recently issued Recommendations of the British X-Ray and Radium Protection Committee 
(October 1948), but it appears that the two conditions of measurement are different in the 
two cases, and over a wide range of radiations the numbers are essentially the same. 
While the recommendations are intended primarily for the protection of the radiation 
worker and not for the patient, the preface very rightly emphasizes the importance of 
caution and restraint in the use of x-rays for diagnosis and the treatment of non-malignant 
conditions. ‘The genetic effects of radiation increase in importance as a larger fraction of 
the population is exposed to ionizing radiation, and diagnostic X-ray examinations become 
more and more frequent. Moreover, the employment of a large number of people in 
industries making use of ionizing radiations calls for, and is being given, serious thought. 
‘The present handbook merely hints at these problems and describes the consultative and 
advisory structure set up in the United States to deal with them. 

Study of the recommendations shows that they are thorough and very detailed. A 
particularly useful appendix gives radiation attenuation curves in lead and concrete for 
beams of x-rays produced at voltages up to 2,000 kv. with worked examples illustrating 
practical calculations of protective requirements. 

The handbook should be in use in all medical institutions using high voltage x-rays. 
It also contains much information which is appropriate to industrial organizations making 
use of these radiations. W. V. M. 


542 


CONTENTS: FOR SECTIONS 


PAGE: 
Mr. Suray N. Gupta. Theory of Longitudinal Photons in Quantum Electro- | 
dynamics. : : : , j ; ; ; ; . 68H 
Prot. Vie, Le Pryce. On the Energy of the Alpha-Particles from Radioactive 3: | 
Nuclei : : : : : . : : , J : . 692 
Dr. W. E. Morrirr. Excited Electronic Levels in Conjugated Molecules—IV : 
Symmetrical Cyanine Dyes. ; ; : f ; : b OG 


Mr. J. Ewies and Dr. C. Curry. Resolution and Analysis of Low Temperature 
Luminescent Spectra of Biand Pb Activated Solids of Simple Crystal Structure 708 
Dr. A. A. Sasry. Influence of Dipole—Dipole eae on the Dielectric Constant 


of Dipolar Substances. 716m 
Dr. P. Wricut. The Effect of Occluded Hydrogen on the Electrical Resistance of : 

Palladium. W927) 
Dr.-R. S. TEBBLE, Dr. I. C. SkrpMore and Dr. W. D. Corner. The Barkhausen 

Peirectn. : : 739 
Dr. K. H. STEwarrt. Domain Wall Movement in a ‘Single Crystal , : 3/68 
Prof. A. Rusinowicz. Sommerfeld’s Polynomial Method Simplified . : 766 
Dr. M. W. Feast. Investigation of the Spectrum of the High Voltage Arc in 

Carbon Dioxide: the CO Flame Spectrum. : ; : : - otk 
Letters to the Editor : 

Dr. F. D. S. BUTEMENT. Radioactive’ t®Er : 2 ae 

Mr. E. R. Raz, Mr. J. G. RUTHERGLEN and Mr. R. D. SMITH. Proton 

Capture Radiation of Fluorine- Lon : : 3 eae oi 
Prof. H. Frou.icH. Isotope Effect in Superconductivity ; 778 
Dr. A. G. Gaypon and Dr. H. G. Wotruarp. Excitation of Spectra i in the 
Inner Cones of Flames ; : ; : ; : : - FIs 

Corrigenda . : : : ; : : : : : : : . 780 
Reviews of Books . : : : : : : : : ; - 730 
Contents for Section B : : ; : : : : : : : “ies 
Abstracts for Section B : : : : : : : : : 2788 


ABSTRACTS FOR SEGIION 


Theory of Longitudinal Photons in Quantum Electrodynamics, by Suraj N. Gupta, 


ABSTRACT. 'The radiation field is quantized by introducing four types of photons— 
two transverse, one longitudinal, and one scalar. The scalar photons are treated by using 
an indefinite metric, and it is found necessary to modify the usual supplementary condition 
slightly. The present theory offers a justification for the symmetrical treatment of the four 
components of the electromagnetic potential, recently applied by a number of authors, and 
proves to be very convenient in applications. ‘The results of physical interest, however, 
are the same as obtained from the ordinary formulation. 


On the Energy of the Alpha-Particles from Radioactive Nuclei, by M. H. L. Pryce. 


ABSTRACT. When the energy of «-particles emitted by radioactive nuclei is compared 
with the energy calculated from the usual semi-empirical formula for nuclear binding 
energies, systematic variations ranging from 4:5 Mev. to —1-0 Mev. are found, which 
exhibit in a very clear way the special importance of configurations containing 82 protons 
and/or 126 neutrons. ‘The general features are readily explained in terms of the shell 
model. ‘The analysis can be used to predict unknown «-energies. Calculated a- energies. 
are tabulated for the heavy nuclei (Z >78) and for nuclei in the rare-earth region. 
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Excited Electronic Levels in Conjugated Molecules—IV: Symmetrical Cyanine 
Dyes, by W. E. Morrirr. 


ABSTRACT. Topics in the electronic theory of the colour of symmetrical cyanine dyes 
have been discussed by means of resonance formalism. Our conclusions for the formamid- 
inium ion have been confirmed by use of the molecular orbital method. As a result of this 
analysis, certain modifications in the qualitative symbolism of colour theory have been put 
forward in the hope that the factors governing absorption may be more clearly understood. 


Resolution and Analysis of Low Temperature Luminescent Spectra of Bi and Pb- 
Activated Solids of Simple Crystal Structure, by J. EwLEs and C. Curry. 


ABSTRACT. Experimental arrangements are described by which. clearly resolved 
spectra have been obtained of the luminescence, at the temperature of liquid air, of seven 
different phosphors of simple crystal structure (cubic) and of mixed crystals. The main 
peaks in each case may be arranged in a band series. ‘The vibration frequency in the lower 
state is nearly the same for the same host lattice with different activators. The energy of 
the electronic transition is determined both by the activator and the host crystal and there 
is a progression of decreasing energy in the series of simple crystals MgO, CaO, SrO, CaS, 
and in the order CaO, CaO-SrO with mixed crystals. The most probable transitions are 
from the lowest level of the upper state. ‘There is evidence that self-excitation occurs and 
for the view that the luminescent centres are large. 


Influence of Dipole—Dipole Coupling on the Dielectric Constant of Dipolar Substances, 
by A. A. Sasry. 


ABSTRACT. Kirkwood’s formula for dielectric constant, €, gives it in terms of a corre- 
-lation coefficient which depends on the probability that any two dipolar molecules have a 
certain configuration. 

The main object of this paper is to calculate the configurational probability of the dipole 
orientations; this is dealt with in the first part of the paper by considering two dipoles and 
teplacing the rest by a homogeneous dielectric. ‘The dipolar molecules are regarded as 
point dipoles and the effect of other forces on the orientation of the molecules is neglected. 

In the second part of the paper the calculation of the correlation coefficient appearing in 
Kirkwood’s dielectric constant formula is carried out. ‘This correlation coefficient is 
evaluated for water, and the value of « then obtained from Kirkwood’s formula is in fair 
agreement with the observed value. 


The Effect of Occluded Hydrogen on the Electrical Resistance of Palladium, by 
P. WRIGHT. 


ABSTRACT. Agreement between the results of theoretical and experimental investigations 
of the absorption of hydrogen by palladium indicates that, at least at higher temperatures, 
the hydrogen is contained in solid solution. Apparatus tor observing the changes in electrical 
resistance of palladium during absorption and evolution of hydrogen is described, with 
which resistance-pressure isotherms from 75° to 150° c. were obtained for pressures 
increasing and decreasing between zero and atmospheric. Factors affecting the resistance 
of the Pd—H system are discussed and an explanation of the form of the isotherms is given; 
a decrease in resistance which is observed, when hydrogen is added to the alloy in the 
B phase, is attributed to completion of the d shells of palladium. 


The Barkhausen Effect, by R. S. Tessie, I. C. SkipMore and W. D. Corner. 


ABSTRACT. An account is given of further experiments on the Barkhausen effect in 
cylindrical ferromagnetic specimens. 'The form of the decay of induction following a 
Barkhausen discontinuity has been investigated and there is reasonable agreement between 
the theory given and the experimental results. As a result of this investigation it has been 
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possible to satisfy the different conditions which are necessary for the detection of the 

discontinuities so as to obtain (a) accurate reproduction of the form of the discontinuities 

and (b) maximum sensitivity in the counting experiments. The number and size of the 

discontinuities are measured, and an estimate is made of the contribution of the discontinuous 

processes to the total change in magnetization over the relevant part of the hysteresis curve. — 
Discontinuous changes in magnetic moment greater than 0-3 x 10~® E.m.lu. (corresponding 

to a minimum volume of 0:9 10-1 cm? for iron) account for from 47% to 86% of the 
total change in magnetization. This is considered in relation to associated work on revers- 

ible changes, and the nature of a Barkhausen discontinuity is briefly discussed with reference 
to other investigations on domain structure. 


Domain Wall Movement in a Single Crystal, by K. H. STEwWaRT. 


ABSTRACT. A single crystal of silicon—iron was shaped so as to encourage the formation | 
of large ferromagnetic domains. Its magnetization curve was found to contain a large 
Barkhausen jump, indicating reversal of magnetization in a volume of the order of 1/10 cm’. 
The rate of change of magnetization was found to be proportional to the excess of the applied 
magnetic field over a certain critical field, and could be made very slow. It is suggested that 
the rate of change is controlled by the braking effect of eddy currents on the movement of 
domain walls. 


sSommerfeld’s Polynomial Method Simplified, by A. RuBINOWICZ. 


ABSTRACT. From two earlier papers by the author it follows that eigenfunctions 
obtainable by Summerfeld’s polynomial method contain either a Riemann P-function or a 
degenerate form of it when the fundamental interval is either finite or infinite respectively. 
‘Comparing the differental equation for the P-function or its degenerate form with the 
differential equation which we get for P from the original differential equation of the 
eigenvalue problem, we obtain standard formulae for the constants determining the eigen- 
values and eigenfunctions. The calculations necessary to give the complete solution of an 
eigenvalue problem which can be treated by Sommerfeld’s polynomial method, are thus 
reduced to a minimum of diffentiations and simple algebraic operations. It is also 
possible, using appropriately chosen P-functions, to obtain the eigenfunctions directly in 
the required form. The method yields also eigenfunctions for eigenvalues in the continuous 
spectrum. 


Investigation of the Spectrum of the High Voltage Arc in Carbon Dioxide: the CO 
Flame Spectrum, by M. W. Feast. 


ABSTRACT. The spectrum emitted by a high voltage arc in CO, at atmospheric pressure 
has been investigated in the region 2000 a. to 9000 a. It was found to consist, mainly, of a 
continuum and the ‘CO flame bands’. A study of the spectrum at different dispersions 
and on photographic plates of various kinds shows that although the Schumann—Runge 
“O, bands do form a part of the spectrum they are relatively weak compared with the bands 
-attributed to CO,. The excitation of the Schumann—Runge O, system in flames is discussed. 
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